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DESTABILIZATION OF FRONTS IN A CLASS OF BI-STABLE SYSTEMS 



ARJEN DOELMAN*t, DAVID IRON*AND YASUMASA NISHIURA* 

Abstract. In this article, we consider a class of bi-stable reaction-diffusion equations in two components on 
the real line. We assume that the system is singularly perturbed, i.e. that the ratio of the diffusion coefficients 
is (asymptotically) small. This class admits front solutions that are asymptotically close to the (stable) front 
solution of the 'trivial' scalar bi-stable limit system ut = u xx + «(1 — u 2 ). However, in the system these fronts can 
become unstable by varying parameters. This destabilization is either caused by the essential spectrum associated 
to the linearized stability problem, or by an eigenvalue that exists near the essential spectrum. We use the Evans 
function to study the various bifurcation mechanisms and establish an explicit connection between the character of 
the destabilization and the possible appearance of saddle-node bifurcations of heteroclinic orbits in the existence 
problem. 

Key words, pattern formation, bi-stable systems, geometric singular perturbation theory, stability analysis, 
Evans functions 

AMS subject classifications. 35B25, 35B32, 35B35, 35K57, 35P20, 34A26, 34C37 

1. Introduction. The class of bi-stable reaction-diffusion equations we consider in this paper 
is given by 



U t = e 2 U xx + (1 + V-U 2 )U 
rV t = V xx + F([/ 2 ,V», 



2 V .^ (1-1) 



where F(U 2 , V; e) is a smooth function of U 2 , V and e such that F(l, 0; e) = and lim e ^o F(U 2 , V; e) 
exists ; r > is a parameter. Thus, the system is such that the background state (U, V) = (±1, 0) 
is always a solution. We furthermore assume that the ratio of the two diffusion coefficients, e 2 , is 
asymptotically small, thus the problem has a singularly perturbed nature. We consider the system 



on the (unbounded) line, i.e. (U,V) = (U(x,t),V(x,t)) with (x,t) £ R x R+. Note that ([TT]) is 
(by construction) symmetric under 

U^-U. (1.2) 



To motivate the structure of (1.1) we introduce the fast variable 



£=-, (1-3) 



so that (1.1) can be written in its equivalent 'fast' form 

U t = % + (1 + V-U 2 )U 



e^rVt = V# + e 2 F{U 2 ,V-e). (L4) 
Since U(x,t) and V(x,t) are a priori supposed to be bounded on the entire domain M x R + , we 



find in the natural (fast reduced) limit, i.e. e — > in (L4), that V = Vb and that U is a solution 



of the well-studied, scalar (standard) bi-stable or Nagumo equation, 

U t = % + (1 + Vb - U 2 )U. (1.5) 

In this paper we interpret the original system, ( |1.1[ ) or ( |1.4[ ), as a scalar bi-stable Nagumo equation 
(|1.5|) in which the coefficient of the linear term is allowed to evolve by reaction and diffusion on a 
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long, or slow, spatial scale. Note that the (slow) dynamics of the V-component are allowed to be 
completely general, except that it is assumed that the full system conserves the symmetry (|1.2|) and 



the background states U = ±1, at V = 0, of the scalar limit (see also Remark 1.1). A priori, one 
expects that the T^-component of front-like solutions will remain small (0(e)) due to the 'boundary 
conditions' V = at ±oo, so that the effect of the slowly varying V(x, t)-component cannot have 
a significant influence on the (well-understood) dynamics of the scalar Nagumo equation. An 
important motivation of the research in this paper is to find out whether or not this intuition is 
correct. 

We will focus completely on the existence and stability issues associated to the persistence of 



the asymptotically stable stationary front solutions of the bi-stable equation (1.5) with Vq = 0. 
In fact, this paper can also be seen as a first step towards analyzing the dynamics (and possibly 
defe cts) of striped patterns in a class of relatively simple bi-stable reaction-diffusion equations, i.e. 



(|l.l| ) for (U, V) = (U(x,y,t),V(x,y,t)) with (x, y) S R . The methods and techniques developed 
in this paper are sup pose d to carry over to the analysis of the existence and stability of spatially 
periodic solutions of (|l . l| ) and their two-dimensional counterparts (the planar fronts and the stripe 
patterns) . See also section [| 

The problem of the persistence of the stable front solution of the scalar bi-stable equation (1.5) 
is quite subtle, as can be expected in the light of recent results on the stability of pulses in singularly 
perturbed reaction-diffusion equations of the Gray-Scott and Gierer-Meinhardt type (J, ||. Such 
systems can also be written in the form ( |l.4|) , however, the scalar limit systems are mono-stable, 
i.e. in essence of the form Ut = U^ — U + U 2 . The pulses correspond in this (fast reduced) limit to 
the stationary homoclinic solution of — u + u 2 = 0. Thus, one would expect that the pulses of 
the full system cannot be stable, since the stability problem associated to the homoclinic solution 
has an 0(1) unstable eigenvalue. Nevertheless, stable pulses of this type do exist in the Gray-Scott 
and the Gierer-Meinhardt equation |Q, On the other hand, the stability of the pulses in these 
mono-stable equations is strongly related to the freedom one has in these systems to scale the 
magnitude of the pulses, i.e. the amplitude of the stable pulses is asymptotically large in e in these 
mono-stable systems. Such scalings are not possible for the fronts in the bi-stable case, since the 
background states (±1,0) are fixed (and 0(1))- 

In the analysis of the front solutions, we will find that it is natural to decompose F(U 2 , V; e) 
into a component that has a factor of (1 + V — U 2 ) and a rest term G(V; e) that does not depend 
on U 2 . Hence, we write (1.1) as, 



U t = e 2 U xx + (l + V-U 2 )U 
rV t = V xx + (1+V -U 2 )H{U 2 ,V;e)+G{V-e), 



><n,j-l ,,. r:v. ( L6 ) 



with G(0, e) = 0. Note that this decomposition induces no restriction on F(U 2 , V; e) since we have 
assumed that F is smooth. In fact 

G(V) = F(l + V, V) and (1 + V- U 2 )H(U 2 , V) = F(U 2 , V) - F(l + V, V). 

We will find that the quantities fp*(0; e) and H(l, 0; e) have a crucial impact on the structure and 
the dynamics of the front-like solutions. Therefore, we define 

dG 

Gi(e) = g^(0:,e) and H = #(1,0; e); (1.7) 

G\ is the main bifurcation parameter used in this paper. Throughout this paper we assume that 
H(U 2 ,V) is non-degenerate, i.e. that H(l + V,V) is not identically 0, and that r = 0(1) (see 



Remark 4.13) 



In section 2 we will show that as long as G\ < and 0(1), the front solutions of (1.5) with 



Vq = persist in a regular fashion, in the sense that the system (1.1) has front solution with 
[/-components that are asymptotically and uniformly close to a front in ( |1.5|) with Vq — 0, and 



with V^-components that are asymptotically and uniformly small (Theorem 2.1). However, if G\ 
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becomes 0(e 2 ) these fronts becom e tru ly singular, in the sense that V becomes 0(1), while the 
[/-component is close to a front of fll.5| ) with Vq ^ 0, on the fast spatial scale (and it converges to 
U = ±1 on the slow spatial scale). Moreover, the front solutions are no longer uniquely determined, 
there can be several types of heteroclinic front solutions if G\ = 0(e 2 ) that may or may not merge 
in saddle- node bifurcations of heteroclinic orbits when G\ is varied (Theorems [2.3| and 2.5). It 
should be noted here that we for simplicity consider G(V) = —e 2 ^V in (1.1) in the singular 
limit G\ — 0(e 2 ) althrough this paper - see Remark We refer to Figure 1 for a numerical 
representation of a regular front (Figure la) and a singular front (Figure lb). The magnitude of 



(a) Regular front, Gi = —1.0. 



(b) Singular front, Gi = -2e 2 



Fig. 1.1. Two asymptoticallly stable front solutions of ( 1. 1 )/( l.t ) plotted on the slow spatial scale x (by a 
numerical simulation). Here H(U 2 ,V) = HoU 2 , G(V) = Gi V , s = 0.1 and Hq = 1. The solid curves represent 
the U -coordinates, the dotted curves the V -coordinates. 

Gi is also extremely relevant in the stability analysis. It can be shown that the (regular) front 
solutions are asymptotically stable as long as G\ < and 0(1) and H + G\ — 2t < and 0(1) 



Theorem 4.3. It seems, at leading order, that the destabilization of the front is caused by the 



essential spectrum, er css , associated to the stability of the front (a css reaches the imaginary axis 



exactly at G\ = or at Hq + G\ — 2t = - Lemma 3.2). However, the analysis also shows that 
there can be eigenvalues near the 'tips' of <r ess , and that it is possible that the destabilization is 
caused by such an eigenvalue, i.e. by an element of the discrete spectrum and not by <7 C ss- These 
'new' eigenvalues do not have counterparts in the (scalar) fast reduced limit problem, they have a 
singular slow-fast nature and may appear through edge bifurcations from the essential spectrum. 

In section 4 we study in detail the nature of the destabilization as G\ < increases towards 0. 
In this section it becomes clear that there is an intimate relation between the geometrical character 
of the singularly perturbed existence problem and the character of the destabilization of the front. 
This is a natural and frequently encountered relation - see for instance jl3| and the references 
there. We establish that a front solution destabilizes at a critical value of G\ — — £ 2 7doubic < by 
an eigenvalue if and only if it merges with another front solution in a saddle-node bifurcation of 
heteroclinic orbits. Moreover, we are able to determine the explicit value of this bifurcation value 
7doubio > 0. If the front does not 'encounter' such a saddle-node as G\ increase to 0, the front will 



be destabilized by a QSS at G± = - see Theorems 4.6 and 4.1(| 



Another way to motivate the analysis of this paper is as follows. In this paper we show that 
the technique of decomposing the Evans function associated to the stability of a 'localized struc- 
ture' (a (traveling) pulse or front) into the product of an analytic 'fast' and a meromorphic'slow' 
transmission function (M ^) can be extended to a class of bi-stable equations. We show that the 
slow transmission function (t2(A, e)) is a natural tool for analyzing the existence or appearance of 
eigenvalues near or from the essential spectrum, and that such eigenvalues play a crucial role in 
the stability of the front. Note that in this sense, the theme of this paper is similar to that of |T4[ ], 
where Evans function techniques are developed to study eigenvalues near a ess in a class of nearly 
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integrable systems. 

The paper is organized as follows. The existence problem is studied in section 2. In section 
3 the basic properties of the linearized stability problem are studied and (the decomposition of) 
the Evans function is introduced. Section 4 is the main section of the paper, in it we develop an 
approach by which the (possible) location and existence of 'slow-fast eigenvalues' near the essential 
spectrum can be studied. This section is split into three parts: a subsection on the regular problem, 
a subsection in which we study an explicit example (G(V) = —£ 2 "f, H(U 2 ,V) — HqII 2 ) in full 
detail, and a subsection in which we study the 'fate' of the regular front as G\ approaches in the 
general case. In section 5 we present simulations which clearly exhibit the impact of the distinction 
between a destabilization by the discrete or by the essential spectrum. Moreover, we discuss some 
related issues and topics of future research. 

Remark 1.1. Large parts of the theory developed in this paper can be generalized to systems 
of the type (1.1)/ (1.4) in which the fast reduced limit system is of the type U t = + B(U 2 ; Vq)U 
for some function B, i.e. to bi-stable systems of a more general nature. We focused on he standard 
case, i.e. B = 1 + Vp — U 2 , since the analysis is more transparent. If one drops the condition 
on the symmetry (|l.2| ), the fronts will in general travel with a certain (nonzero) speed. Although 
the symmetry is used throughout this paper, there is no reason to expect that such asymmetric 
systems cannot studied along the lines of the methods presented here. 

2. The existence problem. We analyze the existence of stationary one-dimensional patterns 
through geometric singular perturbation theory |ll| using the methods developed in || ||. 
Therefore, we write the ODE associated to (1.6) as a dynamical system in R 4 , 



u 


= p 




p 


= -(1 + v 


— u 2 )u 


V 


= eq 




<1 


= e[-(l + 


v — u 2 )H{u 2 , v; e 



(2.1) 



where ' denotes the derivative with respect to the spatial variable £ (1.3) (i.e. £ 'plays the role of 
time'). Note that this system inherits two symmetries of (1.6) 



<,P 



-p,q 



-q and u 



-u,p 



(2.2) 



We consider the 'super-slow' case in which G±(e) — 0(e 2 ) separately in sections 2.2 and 2.3. Note 
that in the fast reduced limit, i.e. e — > in ( [2.1[ ), the monotonically increasing heteroclinic front 
solution is given by (uo,Po, vq, qo), where 



(uo{£,;v ),Po(t;v )) = ( vTTuotanh ( J - + V ° £; 



1 + V 2 

^^sech 

V2 




(2.3) 



and vq and go are constants. 

2.1. The regular case. The main result of this section is, 



Theorem 2.1. Let G\(e) (1.7) be 0(1) and negative. Then, for e > small enough, system 
(2.1) has a symmetric pair of heteroclinic orbits: T^(£;e) = (tt/,(£; s),Ph{£,', e), %(£; e), Qh(£', £)) 
andr~(£;e) = {-u h (£; e), e), v h (£; e), q h (£; e)), with 

lim^ioo T^(£;e) = (±1, 0,0, 0) and lim^ ±co T^(£; e) = (=Fl,0, 0,0); u h (£,; e) and %(£;e) are odd 
and monotonic as functions o/£, Vh(£',s) andphi&s) even. Moreover, |u/,(£; e) — uo(£; 0)| = 0(e) 
(2.S) uniformly on R, ^(^je)!, k/i(£;£)| = 0(s) uniformly onR, and Vh(0; s) is the extremal value 
ofv h (£;e), with 



Vh(0;e) = 



2^GW) 



(l-«g«;0)) H(u 2 ^O),O)d^ + 0(e 2 ). 



(2.4) 
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The orbits (£; e) correspond to the (stationary) front patterns (±t//j(£; e), Vft(£; e)) of (l.t) with 
Uh{£,\s) = Uh(£;e) odd as function of £, V/,(^;e) = = 0(e) euew, lim 5 ^ ±oc £4(£;e) = ±1, 

awe! lim^ ±00 V),(£; e) = 0. 



Proof. As the system is singularly perturbed, we also consider (|2.1[) with the slow scaling 



£ = e£, ( |2.1[ ) is given by, 

evl = p 



ep' = — (l + v — u 2 )u 
v 1 = q 

q' = [-(l + « -u 2 )H(u 2 ,v;e) ~G(v;e)] 



(2.5) 



where ' refers to differentiation with respect x. System ( |2.5| ) is referred to as the slow system. We 
begin by finding the locally invariant manifolds of (|2.5|) in the limit e — ► 0. In this limit, the first 



two equations of (2.5) will reduce to, 

p = Q, -{1 + v -u 2 )u = 0. (2.6) 

The manifold given by (u,p, v, q) — (0, 0, v, q) is not normally hyperbolic and will not be considered. 
However, the manifolds, denoted A4q, determined by (u,p,v,q) = (±vT+~u, 0, v, q) are normally 
hyperbolic and thus by || [yj, ( |2.5[ ) possesses locally invariant manifolds Aif, which are 0(e) close 
toMf. We now determine the leading order correction to theses manifolds. Let the manifold A4f, 
be given by, 



Mf = {u = ±Vl + v + eLP^v, q; e),p = eP ± (v, q; e),v, q} . (2.7) 

To obtain successive approximations of Mf, we can expand U ± = uf + euf + ■ ■■ , and P ± = 
pf + epf + ■ ■ ■ ■ Using the first two lines of (2.5) we find, 



± Q ± dU l 9U X nt \ ± n ± -r- I 2 -r- G ( V '> £ ) tOQ\ 

Pl =27TT^' P2 = ftT^^ G(v;e) ' Ul =0 ' U2 ^ T 4(i + ^/ 2 T (TT^- (2 ' 8) 

Hence, the (slow) flow on the slow manifold is given by, 

v" = -G{v;e) + 0(e 2 ). (2.9) 

To leading order, this flow is integrable. The point (v, q) = (0, 0), that corresponds to (±1, 0, 0, 0), 
is a critical point on A4f. Since G\ < 0, (0, 0) is a saddle on M. f with stable direction (1, >/—Gi) 
and unstable direction (—1, y/— G\). 

A heteroclinic orbit Tf from 0, 0, 0) to (±1, 0, 0, 0) is both an element of W U (MJ) and of 
W s (A4f). Here we will only consider Tf. The existence of follows from the symmetry (|2.2| ). 
The orbit remains exponentially close to W u (— 1, 0, 0, 0)]^- before it 'takes off' and a makes 
a 'jump' through the fast field. After that, it 'touches down' on A4f and remains exponentially 



close to it (and to W U (1,0,Q,0)\ M + - see Figure |2.l| . The change in q by the passage through 
the fast field is 0(e) fl2.l| ), therefore must take off from M~ and touch down on A4f with a 
g-coordinate that is 0(e). Since is asymptotic to the saddle points (0, 0) S M. f, it follows that 
the v-coordinate of must also be 0(e). Note that we have used here implicitly that G\ — 0(1). 

We will determine whether such a trajectory, as T~f, is possible using a Melnikov method. Both 
W U (M~) and W s (Mf) are 0(e) close to the family of heteroclinic orbits in the fast reduced limit of 



(2.1) given in (2.3). The leading order distance between W U (M~) and W s (A4f) can be determined 
by a Melnikov function for slowly varying systems JLsfl . Both W U (M~) and W s (Mf) intersect 
the hyperplane {u — 0} transversally. Note that W s,u (Mf) H {u = 0} is 2-dimensional, thus, since 
{u = 0} is 3-dimensional, one expects a 1-dimensional intersection W u (M.~) <1W S (Mf) D {u = 0}. 
The separation between W U (M~) and W s (Mf) is, at leading order, measured by the integral, 



6 



A. DOELMAN, D. IRON AND Y. NISHIURA 



Here the wedge product refers to the scalar cross product and g§ solves the differential equation, 
<7o£, 

the leading order splitting distance, 



s| (if) = ' §K^) = ^' Substituting (2.3) into (2.1C) results in the following expression for 



A = - 



sL 

V2 



tanh 



V2 



sech 



V2 



d£ =-q V2. 



Thus, W U (M~) n W s (Mj) n {u = 0} must be 0(e) close to q = 0. By the symmetries ( |2.2|), w e 
conclude that W U {M^) H W(.M + ) n {u = 0} must be identically q = 0. Hence, again by (§J), 
any solution that connects A4j to .M^~ must have a u component that is odd with respect to £ 
and a i> component that is even with respect to £. 

We are now ready to determine the take off (touch down) curves T~ C (T£ C M.f) || [| . 
These curves represent the points at which the one-dimensional family of orbits in W U (A4~) n 
W 8 {M^) leave (land on) Aif . Let the elements of this family be denoted 7(£;p), where the 
parameter p > corresponds to the p-component of 7(£;p) as it crosses through {u = q = 0}. 
Note that the 7- family forms the Fenichel fibering of W U (M~) f~1 M /s (tW+) |J and that each 



7(£; p) is asymptotically close to an unperturbed orbit given in ( |2.3|) . To each 7(^5 p) we associate 
two orbits, j M -(&p) C .Mj and 7x+(^;P) C .M+ by the fact that \\^;p) - J M ±(&p)\\ is 
exponentially small if ±£ > 0(e _1 ). We define T~ and Ti" as the collections of base points of the 
Fenichel fibers on M.~ and on Mz , 



T o = (J ^M-i^P) > 



U 7m+(°;p)- 



(2.11) 



p>0 



p>0 



We can compute the leading order structure of T~ and Tf by considering the effect of the journey 
through the fast field on the slow variables v and q. Since = eq and q = 0(e) it follows that 
the change in v through the fast field is of higher order, i.e. 0(e 2 ). By construction, q will be an 
odd function of £, thus the value of q for a given v on T~ must be — ^Aq(v), where Aq(v) is the 
change in q due to one full pass through the fast field (during which v remains (at leading order) 
constant, v = t>o). Similarly, the value of q on must be ^Aq(v). Si nce we alr eady know that 
both v and q must be 0(e) in this regular case, we compute Ag(0) (by (2.1), ( [2.3| )), 



/OO />C 
$i„=ode, =-e / 
-OO •/ — ! 



1 - tanh" ( 4= 
^2 



//Haulr ( X ) .(„,/(, 



0(e 2 ). 



To establish the existence of the heteroclinic orbit r^(£), we consider the intersection T Q n 



W^(-l, 0,0,0)| 
given by {g = 
position of T~ 



m: 



on M\ 



0(e) close to (-1,0,0,0) 
0(e 2 )} and {g 



Thus, T~ and W u (-1,0, 0,0)1^ are 
V^G~iv + 0(e 2 )}. 



with ^"(-1,0,0,0)1^- and of T 



T~ n W u (-1, 0,0,0)].,- is given in 



(0) 



Figure |2.1| shows the super- 
with 1^(1,0,0,0)1^+. The w-coordinate of 



Thus, we have established the existence of an orbit 
T+ g W U {M~) n W s {Mt) that is asymptotic to (-1,0,0,0) G M~. Since T+ passes through 



{u = 0, g = 0} during its jump through the fast field, it follows by the symmetries (2.2) that 
is indeed the orbit described in the statement of the Theorem. As was already mentioned, the 

□ 



existence of T h also follows immediately from ( |2.2| ). 



Remark 2.2. We note that if G\ = 0(e CT ) for some a g [0,2) then the intersection of T a 
and 1,0,0, 0)1^- will result in a value of v Q of 0(e 2 -< T ) < 0(1) Q, thus r+(£) will still 

be a regular perturbation of the orbit in the scalar limit. Moreover, this argument also shows that 
singular orbits may exist for G\ — 0(e 2 ). 

2.2. The super-slow limit: an example. In this section we consider the 'significant de- 
generation' Gi(e) = 0(e 2 ). For simplicity, we only consider the case in which the flow on the slow 
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W u (-1, o,o,o)\ Mi 




M: 




^(1,0,0,0)|^ 



m: 



Fig. 2.1. Superposition of the take off and touch down curves T Q d with W s ' u (±l, 0, 0, OJI^i ■ The intersections 
T~ nW u (-1,0,0, 0)\ M - and T+ n W s (1, 0, 0, 0)| M + determine the heteroclinic front solution e). The dotted 

arrows inidicates the orbit (§) 'taking off ' and 'touching down'. 



manifolds Aif is linear, i.e. G(v;e) — e 2 Gi(e) == — £ 2 7, where 7 does not depend on e. Moreover, 
we first consider an explicit expression for H(u 2 ,v; e), H{u 2 ,v;e) = Hpu 2 . T he case of a general 
H(U 2 , V) will be considered in the next subsection. We refer to Remark 2.4 for a brief discussion 



of the case of a general function G{V). System ( |2.l| ) reduces to 




U > (2.12) 
— u 2 )H u 2 + e 2 jv] . 

This system has various types of (singular) heteroclinic orbits. 

Theorem 2.3. Assume that G(V) = —e 2 ~fV, H(U 2 ,V) = H U 2 and that e is small enough. 



(i): H a > 0. If 7 > 7doubic, where 7doubic = %Hq + has two pairs of hetero 



clinic orbits, T~^' J '(£; e) = (w^ (£ ) , (£) , v^(£), g^(£)); 3 — 1)2, and their symmetrical counterparts 
r^e) = (-«l(0 ? -jteX(0»9h(0)> withlim^ ±oo r+' j (Z;e) = (±1,0,0,0). In the fast field 
u h(0> respectively v 3 h ({;), is asymptotically and uniformly close to Ua(£\Vj) (pT3j resp. Vj,- the con- 
stants Vj are the zeros of s/l v = ^V2Ho(v + l) 3 / 2 so that < v\ < 2 < V2 (at leading order). In 
the slow field, T^' J (£;e) is exponentially close to W u,s (±l, 0,0, 0)\ M ± C Mf . The orbits T^' 1 (£; e) 
and T h ' 2 (£;e) merge in a saddle-node bifurcation of heteroclinic orbits as 7 J, 7doubie- There are 
no heteroclinic orbits for 7 < 7doubic- 

(ii): H Q < 0. The relation ^jv = jV2H (v + l) 3 / 2 has a unique zero for all 7 > and there 
is one pair of heteroclinic orbits r^(£;e) for all 7 > 0. These orbits have the same structure as 
described in (i). 



The orbits F^' correspond to the front solutions (Z7^'"^(£; e), V^'^ (£; e)) of \l.d\} with 

Ufr (£>e) — ± u h(^! £ ) and Vj^^'fce) = v k (^;e) even as function of £. 

Proof. The essence of the analysis of the super-slow system is similar to that of the regular 
case. The important difference being that, although the change in q by a 'jump' through the fast 
field is still 0(e), the ^-coordinate of the heteroclinic orbit may now be 0(1), due to the super-slow 
character of the flow on J\4f. It is this difference that will cause the bifurcation and the formation 
of the second orbit in case (i). The flow on the slow manifold is now 0{e 2 ), i.e. super-slow, and is 
at leading order governed by, 

v" = e 2 jv. (2.13) 



Since the right hand side of this equation is 0(e 2 ), one might expect that one needs to incorporate 
the higher order corrections to the approximation of J\4f (2.8) to determine the leading order flow 
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on A if. However, the 0(e 2 ) correction contains a term with a q 2 factor and a term with G(v) 
(2.8). Since we consider q — 0(e) on A4f and since G(v) — 0(e 2 ), the resulting correction will 



not be of leading order. 

Again the equilibria on Aif are saddles, with stable and unstable directions, (±1,6^/7). As 
in Theorem [2.l| we only consider the orbit that jumps from A4~ to A4f (the others follows from 
the symmetry (2.2). We repeat the Melnikov calculations and again conclude that, W U (A4~) n 
W s (Mf) n {u — 0} must be identically q = 0. Hence, again by (2.2), any solution that connects 
M.~ to M. J must have a u component that is odd with respect to £ and a v component that is 
even with respect to £. 



Wc define the take off, T~ , and touch down, Tt, curves as in (2.11). We find the leading 
order behavior of and T~ , by calculating the change in q as we traverse the fast field. As 
in the regular case, v remains a constant up to 0(e 2 ) and the value of q on the take off (touch 
down) curve must be — ^Aq(vo) (^Aq(vo)), where vo is the (leading order) constant value of the 
u-coordinate of the orbit that is heteroclinic to A4f in the fast field. The calculation of the c han ge 
in q is similar to that of the regular case except that vq now effects the leading order term (2.3), 



Aq(v a ) = - £ i/ (l + ^) 2 /r oo [l-tanh 2 ( v ^±i£)]tanh 2 ( 
- _^ iJo(wo + 1) 3/ 2 + 0(£ 2 ) 



^±±£ d£ + 0(e 2 ) 



= —e 



The heteroclinic orbits are again determined by T Q n W u (— 1, 0, 0, 0)\ M -, where T Q = {q = 



M- M- 

(a) (b) 

Fig. 2.2. Superposition of TJ~ with W u (— 1, 0, 0, 0)| in the super-slow case with Hq > 0, for 7 > 7doublc 
(a) and 7 > 7 doub i c (b). 

~±Aq(v ) + 0(e 2 )} and W u (-1 7 0, 0,0)\ M - = {q = e^v + 0(e 2 )} 

^o(«o + l) 3/2 -V7«o, (2.14) 
see Figure E^3. Thus, in the super-slow case, a heteroclinic orbit may leave A4~ with a w-coordinate 



of 0(1). Now if H > and 7 > 7doubic = I-Hq + 0( £ )j (2-14) has two possible solutions, vq = 



v 



3' 



j = 1, 2, with < vi < 2 < t>2 (at leading order). These intersection correspond to the heteroclinic 
orbits r^ J (£). For 7 < 7doubic, there are no solutions to (2.14) and thus no heteroclinic connec- 



tions exist: the orbits T^ ,1 (£) and r^"' 2 (£) have coalesced at 7 = 7dout>ic- In the case that H < 0, 



(2.14) has a unique solution for all values of 7 > 0, there is only one pair of heteroclinic orbits. □ 



Remark 2.4. If G(V) is not linear in the singular limit (i.e. G\ = 0(e 2 )), then the analysis 
becomes more involved, but there are no essentially new phenomena. In this case, the magnitude 
(w.r.t. e) of the second derivative of G(v) at v — will start to play a role comparable to G\. 
Moreover, the flow on M~ is nonlinear, so that VP"' S (±1, 0, 0, 0)1^- is no longer a straight line 
(at leading order), therefore, many 'new' intersections of T~ n W v (—1, 0, 0, 0)1^- , and thus 'new' 
heteroclinic orbits, may appear. 
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2.3. The sup er-slow limit: the general case. Wc now consider the general su per-sl ow 
problem, i.e. (2.1) with G = — s 2 -fv. The treatment of the general super-slow case and (2.12) is 
in essence identical to that of the previous section. However, the statement of the main results 



cannot be formulated as explicit as in Theorem 2.3, as long as there is no explicit expression given 
for H(U 2 , V). Nevertheless, the character of the existence result is similar to that of Theorem 2.3, 
there can be various kinds of heteroclinic orbits that might coalesce in saddle-node bifurcations. 



As in the proofs of Theorems 2.1 and 2.3, the existence of the heteroclinic orbits is established 



by the intersection of T a and W u (— 1,0,0, 0)\ M - , i.e. by the solution vq of 



= 77 



[l + v - Uo(£;u )] H(ul(€;vo), v )d£, 



(2.15) 



at leading order. Note that the right hand side equals — ^Aq(vo), i.e. half the accumulated change 
in q during a circuit through the fast field, and that we have used ([2~3|). 

Theorem 2.5. Assume that G(V) = —e 2r fV, and that e is small enough. System ([^ 
has n > pairs of heteroclinic orbits, r^ ,:, (^;£) = (iu^' 3 '(£), ip^' 3 '(£), v h ' 3 '(£), q^ J (£)), where 



J 



with lim 



T h J (£; e) = (±1, 0, 0, 0). The number n — n("f) is given by the number 
of solutions Vj of (2. It). In the fast field u 3 h (^), respectively * s asymptotically and uniformly 

close to uo(^;vj) (2.S) resp. Vj, where the constant Vj is the j-th zero of (2.15). In the slow field, 
is exponentially close to W u ' s (±1,0,0, 0)\ M ± C 



1 h 



Two orbits T h ' J (^;e) and T h '■ 7+1 (£; e) coalesce in a saddle-node bifurcation of heteroclinic orbits 



at a certain value 7 = 7d oub i c , */ ^ e zeroes Vj 



< v 



•3+1 



of ( 2.11 ) merge, i.e. if the intersection 



r-nw(-i,o ) o,o)| A4 - 



-transversal. 



The orbits rJ J (£;e) correspond to the front solutions (U* ' 3 (£; e), J (£; e)) of ( \l.6\ ) with 



— ±u° h (^;e) odd, and V^ ,3 (^;e) — v 3 h (^;e) even as function of £ 



r±,j 



The proof of this result is in essence identical to that of Theorem 2.3. In Figure 2.3 



two 



examples of the possible richness of the intersection T Q n W u (— 1, 0, 0, 0)| w- are given 




Fig. 2.3. Two examples of the possible character of the intersection T n W u {— 1,0,0, 0)\ M - for a given 
H(U 2 ,V); (a) there are 3 different singular heteroclinic orbits (b) 4 heteroclinic orbits. 



3. The stability of fronts. With a slight abuse of notation we (re-)introduce u(£) and 



by 



U(t, t) = U h (C, e) + u(t:)e M , V({, t) = V h (£, e) + 



substitute this into ( jl.6| ), and linearize 

u Si + [1 + V h - 3C/2 _ x)u 



-U h v 



{2 



H{Ul V h ) -(1 + V h - Ul)j&(Ul, V h )} U h u 

H{Ul, Vh) (1 + V h Ul)^§{Ul,V h ) + ^(V h ) - rX] v} 



(3.1) 
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Note that the front pattern (Uh ( Q, Vh(Q ) corresponds to any of the regular or singular heteroclinic 



bits r, '•'(£) of Theorems 2.1, |2.3|, 2.5. In the stability analysis of forthcoming sections we will 



only consider the front patterns of +-type, i.e. those fronts for which lim^-too t//j(£; e) = ±1. 
Thus, we do not explicitly consider their symmetric counterparts. Due to the symmetry (1.2) this 
is of course also not necessary. The coupled system of second order equations (3.1) is equivalent 
to a linear system in C 4 , 

C =A(fcA,e)0 with 0(0 = (u(0,J>(0> "(0.9(0), (3-2) 
where A(^; A, e) is a 4 x 4 matrix with Tr(A(£; A, e) = 0, and = p, vg. = eq. It follows that 



lim A(t;\,e) d ^ A±(\e) 

£— >±oo 



/ 1 \ 

2 + A Tl 

e 

\ ±2eH a -e(H + G x - rA) / 



(3-3) 



(1.7). The matrices ^4^, have the same set of eigenvalues Aj(A, e), i — 1, 2, 3, 4, 
A?, 4 (A, e) = A + 2 + 0(e 2 ), A| 3 (A, e) = £ ^ 2 - A( Gl + ^- 2r) - 2 Gl 



+ 0(e 4 ). (3.4) 



Note that both expansions break down as A approaches —2 (see Remark 3.1). We define a branch 
cut such that for z £ C arg(y / z) £ (— iw, hn], so that the Aj's can be ordered 



Re(A 4 (A,e)) < Re(A 3 (A,e)) < < Re(A 2 (A, e)) < Re(A 3 (A,e)) 



(3.5) 



This ordering of course breaks down if A £ (T css , the essential spectrum associated to ( |3.1[ ) / ( pT2|) , 
since a css coincides with values of A for which either Re(Ai i 4(A, e)) = or Re(A2 i 3(A, e)) = M; 
see also section 3.1). The eigenvectors E^(e, A) of the matrices A^(e, A) associated to Aj(A, e) are 
given by 



1 

Ai, 4 (A,e) 
0(e 2 ) 

2H 



/ ± J L- 2 +0(e 2 ) \ 
0(e 2 ) 
1 

V iA 2 , 3 (A, e ) J 



(3.6) 



(for A + 2 = 0(1) - Remark fy). 

Remark 3.1. The expansions ( |3.4| ) and (3J3) are only valid for A + 2 = 0(1) with respect 
to e. It is straightforward to check that A 2 4 (A, e) — 0(e) = A| 3 (A, e) if A + 2 = 0(e) and that 
in general, when A + 2 = 0(e CT ) for some a £ [0, 1], A? 4 (A,e) = 0(e CT ) and A 2 . 3 (A,e) = 0(e 2 ~ CT ). 
Thus, Ai j4 cannot be assumed to be large/fast compared to A 2 ,3 if A + 2 = 0(e). Since A » — 2 
is way into the stable region, we do not consider this degeneration further and assume throughout 
this paper that |Aa 3 j <C |A 14 |. 

3.1. The essential spectrum. The essential spectrum associated to the stability of the front 
patterns (U, V) = (Uh(£), Vh (0) * s ^ uu y determined by the spectrum of the linear stability problem 
for the (trivial) background states (at ±00) (U, V) = (±1,0) [0. Therefore, we introduce HI 
and a, f3, A £ C by 



U(x,t) = ±1 + ae lki+xt , V(x,t)=0> 



and substitute this expression into (1.6) (using (1.3)). This yields the matrix equation 



f-k 2 ~2 ±1 
\T2e 2 H -k 2 + e 2 {H + G x )) K/3 



= A 



a 
e 2 t/3 
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where G\ and Ho have been introduced in ([l/i]). Thus, A = A(fc 2 ) is a solution of the characteristic 
equation 

Q(X, k) = (A + k 2 + 2)(e 2 rA + fc 2 - £ 2 (FT + d)) + 2e 2 H a = 0. (3.7) 



Note that this equation holds for both background states (±1,0), due to the symmetry (L2). We 
may conclude 

Lemma 3.2. The essential spectrum cr css associated to \3. i[ ) is given by the solutions A = Xik 2 ) 
of §3~7\ ) with k e R; <r css is stable, i.e. a CSB € {Re(A) < 0}, if Gi < and H + Gi - 2t < 0. 
Proof. The two conditions in this lemma are obtained directly from 



Ai + A 2 = ± [e 2 (H + G 1 -2 T )-k 2 (l+e 2 T)] < Vfc 
AiA 2 = £ [k 4 + k 2 (2-e 2 (H + G 1 ))-2e 2 G 1 ] > V k. 



(3.8) 



Both relations attain their extremal value at k — 0. □ 
However, we need to have more information on the essential spectrum than just this stability result. 
In section || we will see that the appearance of edge bifurcations is closely related to the structure 
of (les s- W e focus on the stable case Gi < and Ho + G\ — 2r < 0. It is straightforward to check 
that (3.7) has two solution Ai^A;) € K for all k S K if Ho < 0. As Ho passes through zero two 
fc-intervals, (— k + , —k~) and (k~ , fc + ) (0 < k~ < k + ) appear in which Ai j2 (fc) are complex valued. 



These intervals merge (i.e. k { 0) as H approaches (v2t — yJ—G±) . For (v2r — y/—Gi)~ < 



H < 2t — Gi (which is a non-empty region), Ai j2 (fc) G C if — fc+ < k < k + . See Figure 3.1 



3.2. The Evans function. The use of the Evans function in the analysis of linear systems 
associated to the stability of traveling waves is by now well-established. Here, we give a brief 
exposition of the characteristics of the Evans function in reaction-diffusion systems. We refer to 
|l], [It], [l(], Q U for the full analytic details of the statements in this section. 

We define the complement of the essential spectrum by 

C e = C\a css . (3.9) 

For A € C e the ordering ([T|) holds, so that 

Lemma 3.3. For all A € C e there exist two two-dimensional families of solutions $_(^;A,e) 
and $+(£; A,e) to $J\ ) such that lim^ ±00 0±(£; X,e) = (0, 0,0, 0)* for all 0±(f; X,e) S $±(£; A,e); 
^±(£l A, e) depend analytically on A. 

An eigenfunction of ( |3.2[ ) must be in the intersection of <!>_(£; A, e) and <£+(£; A, e). Since 
Tr(A) = 0, we therefore define the Evans function T>(X, e) by 

V(X, s) = det[Md; A, e),M(; A, e), <M£; A, e), <M£; A, e)], (3.10) 

where {0i,</> 2 } (respectively {4>3,4>4}) span the space $_ (resp. $+)■ The Evans function is 
analytic in A e C e , its zeroes correspond to eigenvalues of ( |3.2| ) counting multiplicities |l|, |l7j ]. 
Of course, this definition does not determine £>(A) uniquely. However, this can be achieved by 
choosing </>i(£) and 2 (£) as follows 

Lemma 3.4. For A £ C e i/iere is a unique solution (f>i(£;X,e) € <!>_(£; A, e) o/ (S.lj smc/i 

lim </»i(^;A,e)e- Al(A ' £)? = £f (A,e). 

4 — > — oo 



(3.4), (3.6). There exists an analytic transmission function tx(X,e) such that 

lim < ? !) 1 (^A,e)e" Al(A ' e) « = *i(A, s)E+(X, e). 
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For X € C e such that ti(X,s) there is a unique solution 4> 2 {£,] A,e) € | I > _(£; A,e) of (3.i), that 
is independent of <fri(£; A, e), that satisfies 

lim cj) 2 (£;\,E)e~ A2 ( x ^ = E^{X,e) and lim 2 (£; A, e) e - Al{x ^ = (0, 0, 0, 0)*. 

There exists a second meromorphic transmission function ^(A, ff), i/iai is determined by 

lim 2 (C;A,e)e- A2 ( A ^ = t 2 (A,e)^(A,e). 

The solutions 03 4(^5 A, e) £ A,e) of (3.2) can be denned likewise. Since Y^t=x^i(^ £ ) = ® 

V(X,s)=dct{M0e~ Ali ,M0e- A2i ,U0e- A3i ,M0e- Ai % 

so that T>(X,e) can be decomposed into a product of ii(A, e) and t 2 (X 7 e) by taking the limit 
6^ +00. 

Lemma 3.5. Let X e C e , then 

V(X, e) = ti(X, e)t a (A, e) det [£+(A, e), £+(A, e), E+(X, e), E+(X, e)] . (3.11) 



We conclude that the eigenvalues of ( p.2|) correspond to zeroes of the transmission functions t\ (A, e) 
and ^2(A,e). However, we will see that a zero of ti(X,e) does not necessarily correspond to a zero 
of T>(X,e), since t 2 {X,e) can have poles (see also H). 

3.3. The fast eigenvalues. Next section will be devoted to the analysis of (the zeroes of) 
t 2 (X,e), here we consider the zeroes of the fast transmission function ti(A, s). In order to do so, 
we first consider the stability problem associated to the front solution E//(£; Vb), with Uf(£; Vq) — > 
±(1 + Vq) as £ — > ±oo, of the scalar fast reduced limit problem (|l.5|), 

«>« + (1 + V - 3 M 2 (£; V Q ) - X)w = 0, (3.12) 

since £//(£; Vq) = wo(£; K)) (2-3). This system can be written as a linear system in C 2 , 

^ = B(& A)V with ^(0 = (u(0,p(0). (3-13) 

where £?(£; A) is a 2 x 2 matrix of which the coefficients are by construction 0(e) close (uniformly 
in Qto those of the 2x2 block in the upper left corner of the 4x4 matrix A* (£; A, e) defined 
in (3.2), if we set Vq = 14,(0). The Evans function associated to this prob lem can be written 
as T>f(X) = det[-0i(£, A), ip^fe, A)], in which an d ip4,(Q are solutions of ( |3.2| ) determined by 

lim^-oo j> iX£)£ VxT * i = (1, VX+ly and lim^oo ^(^e^ 3 * = (1, -a/XTI)* (where ±y/X + 2 
and (1, ±y/X + 2)' are the eigenvalues and eigenvectors of the matrix Boo(X) = lim^ioo £?(£; A) 
(compare to ( |3.4[ ), ( |3.6| ))). As for the full system, we can define an analytic fast reduced transmis- 
sion function t f {X) by lim^oo ^i(£)e _v/ *+^ = i/(A)(l, vX+2)*, so that 



V f (X)= lim det[V^ (0.^4(01 =det[t / (A)(l,VAT2) t ,(l,-v / AT2) t ] = -2t f (X)VX + 2. 

£ — ''CO 

The transmission function ti(A) is, by construction, asymptotically close to its fast reduced limit 
t f (X). 

Lemma 3.6. Let X{ £ C e such that tf(X{) = 0. There is a uniquely determined Xi(e) with 
lim e ^o Aj(e) = X{ such that t 1 (X i (e),e) — 0; t 1 (X,e) ^ for X ^ A,(e). 

The proof of this Lemma is completely analogous to the proofs of similar statements in Jj], 
[To[ [|, |). Hence, we find (the leading order behavior of) the zeroes of £i(A, e) by computing the 

spectrum of ( 3.1 2| ) . By ( [2~5| ) and by introducing r\ = \/|(l + Vb) we can write ( |3.12j ) as 



/ 6 9 \ , 2A 

- 2 P 2 )w = with P 2 = — + 4, 

V cosh n J 1 + Vq 
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which is a well-studied problem of Schrodinger/Sturm-Liouville type (see for instance J2(], [|). It 
has discrete spectrum at P — 1 and P = 4 and essential spectrum for P € iR. We conclude that 
the eigenvalues of (3.12), and thus the leading order approximations of the zeroes of ii(A), are 
given by 



A 



0, A^ 



- 5 (l + %)<0. 



(3.14) 



The essential spectrum of (3.12) is given by 

*L = {A < -2(1 + V )}. 



(3.15) 



We conclude this subsection by stating two simple, but useful results: 

Lemma 3.7. Let (it(£; e), ?;(£; e)) be a pair of eigenfunction solutions of (3.1) associated to a 
simple eigenvalue X(s), then either m(£) is even as function oft; and i>(£) odd, or u(£) is odd and 
i>(£) even. 

Proof. We write (3^) in the following way, 



v^=e 2 [F o (0u + F e (0v}. 



(3.16) 



By construction, Uh is an odd function of £ and Vh is an even function of £. It thus follows that 
the above functions, F Q and F e must be odd and even functions of £ respectively. Let (u, v) be an 
eigenfunctions associated to the eigenvalue A. We decompose (u, v) into odd and even components, 
u = u + u e , v = v + v e , where u , v are odd and u e , v e are even. By the parity of the functions 
Uh, Vh, F and F e it is clear that (u ,v e ) and (u e ,v ) form two independent solutions of the 
eigenvalue problem associated to the eigenvalue A. Since we have assumed that A is simple, we 
have a contradiction. □ 



Lemma 3.8. Assume that the eigenfunction solution v(£) of (3.1) with eigenvalue A(e) is odd, 
then A(e) = ; so that (u(0,«(0) = (Uh,d& e), Vfc,e(& 

We will see in section ^ that there can be several eigenvalues for which u(£) is odd and t>(£) 
even. 

Proof. It is clear that there is an eigenvalue A = associated to the derivative of the front 
(u(£),i;(£)) = (Uh^{^', £),Vh^(^; e)). We assume there is another eigenfunction with v odd. Since 
Vff is 0(e 2 ) and v is odd, it follows that |v| <C 1 on the fast spatial scale. Hence, the equation 
for the it-component is to leading order homogeneous and given by ( 3.12|) (with u replaced by w). 
Lemma 3.7 implies that u is even. Since the only even eigenfunction of (3.12) is Uh£ with eigenvalue 
0, it follows that the leading order behavior of u is given by Uh,£ and that A is asymptotically close 
to 0. We thus write, 



U h ,f + S(e)ui 



S(e)vi , A = S(e)X(e) 



(3.17) 



where 5(e) — > as e — > and A(0) ^ (i.e. S(e) represents the leading order magnitude of A). We 



substitute (3.17) into (3.1), to get the following equation for u\, 



*i >€£ + (1 - 3E#)ui = Al7 h)£ - U hVl . 



This equation has the solvability condition, (\Uh,£ — UhV\)Uh,£ dy = 0. Now, v\ and Uh are 

odd while Uh,£ is even, thus we have that A = 0, contradicting our assumption. So the only pos- 
sible eigenfunctions with v odd must correspond to a eigenvalue and hence, (u, v) = (Uh,£, Vh,e). □ 
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4. Slow-fast eigenvalues and edge bifurcations. The 'slow-fast eigenvalues' are the eigen- 
values that exist due to the interaction of the fast [/-equation and the slow U-equation in (1.6), 
thus, these eigenvalues do not have a counterpart in the fast reduced scalar limit problem (1.5). 
The slow-fast eigenvalues correspond to the zeroes of the f2(A;e), since this transmission function 
is based on a balance between slow and fast effects. See also Remark 4.5. 

In order to study the combined effect of slow and fast terms, we need to define the region in 
which the fast £-jump takes place more accurately 



{ £e(--L -L)} or {set-Vi.Vi)} 



(4.1) 



(1.3). Note that the exact choice of the boundaries of If is not relevant, any choice will be suitable 
as long as it is in the transition zone between x and £ (i.e. on the boundary of If we must have 
\x\ < 1 and |f | > 1). 



4.1. The regular case. Again, we first consider the case G\ = 0(1) (1.7). In the slow 
coordinate x (1.3), i.e. outside the region //, the equation for u reads 



(1 - Wl - X + 0(e))u = -U h v + 0{e 2 u xx ) 



(4.2) 



(3.1), since Vh(tj) — (e) on M (Theorem p.l| ). Thus, u can be expressed in t erm s of v outside the 
fast ^-region If (4.1). Using that L^(f;e) = 1 + 0(e) outside If (Theorem 2T), we find for the 
^-equation of (3.1) on the slow a;-scale, 

2H(1, 0)U h + 0(e)] u - [H{1, 0) + |f (0) - rA + 0(e)] v 



2Ho 
A+2 



Ho-d +rA + 0(e) v + 0(e 2 v xx ) 



(1.7). Hence, outside If 



-# A + A(A + 2)t-G 1 (A + 2) 



A + 2 



0(e) 



(4.3) 



uniformly in f . The w-equation is thus at leading order of constant coefficients type. By (3.7) and 
(3^) we have on the f -scale 



Q(A;0) 



A + 2 



0(e 3 ) 



v= [Al 3 (X,e) + 0(e 3 )]v, 



(4.4) 



In order to determine an expression for £2 (A, e), we need to control the solution <p2 (£; A, e) (Lemma 
[T4|) of (3.5). This is done in the following Lemma. 

Lemma 4.1. For all X G C e such that t\{X,e) 7^ there exist 0(1) constants C_, C+ > and 
a third meromorphic transmission function t 3 (X,e) such that 



fo(£;X,e) 



[£ 2 ~(A) + 0(e)] e ^WZ +0(e c -Z) for £ < 
i 2 (A)£+(A)e A2 ( A )« + t 3 (X)E+(X)e A3( - x ^ + 0(e~ c +«) for £ > ^ 



(4.5) 



Moreover, there exists an 0(1) constant Cf such that \\4>2{C)\\ — C/ f or £ € If- The v-coordinate 
of 02(C) satisfies v(£) = 1 + 0(y/e) on If, so that 



ta(A,e)+t 3 (A,e) = 1 + 0(Vsj. 



(4.6) 



Proof. The behavior of <p2(0 outside If is determined by (4.4) and ( |4. 2| ) . The approximation 
(4.5) for £ < —1/^/e follows from the definition of 02(C) (Lemma 3.4). This same Lemma establishes 
the leading order term in (4.5) for £ — > 00. The transmission function ts(A, e) measures the 
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component of 2 (£) tha t decays on t he slow spatial scale x. In side //, — 0(e 2 ) (3.1) and 
AjjjAjE) = 0(e 2 ) (3.4), so that (4.6) follows. As in section 3.2 we refrain from giving the full 
analytic details of this result, since these are essentially the same as in □ 
The transmission function i 2 (A,e) can be determined by the methods originally developed in ||. 
We deduce from Lemma 4.1 that the total change in over // is given by 



A slow « ? = 2e(t 2 (A) - 1) 



Q(A;0) 
A + 2 



(4.7) 



This change in must be an effect of the evolution on the fast £-scale, that is given by 



Af ast l>£ 



"«|{«=«, n ,^=i}^ + 0(e 2 Ve), 



where Ui n (£) is a bounded solution of the inhomogeneous problem 

u« + (1 - 3C/ 2 (£; 0) - X)u = -U h (t; 0) 



(4.8) 



(4.9) 



(rec all tha t v(£ ) = l + 0(y / e) in//). The transmission function i 2 (A; e) is determined by combining 
Since, a priori A s \ ow v^ = 0(e) and Af ast V£ = 0(e-Je) we are led to the following 



(|4jj) and (|4 
conclusion. 



Lemma 4.2. Consider A <E C e n{Re(A) > — 2+6} for some 5 > independent ofe. Let A 



/ 



'it the second eigenvalue of the limit system ) with Vq — (3.14), and let A + (0) and A (0) be 



(3.7). Then, t 2 [X) =l + 0(^i) if |A - \{ |, |A - A+(0)|, |A - A"(0)| 



) if\\-\ f 2 \ = 0(e°), |A-A+(0)| =0(e 2CT ), or |A - A~(0)| = 0(e 2ff ) for 



the solutions of Q{\, 0) 
0(l);ta(A) = l + 0( £ 2- 
some g £ (0, \). 

Thus, this Lemma establishes that i2(A,e) can only be zero in {Re(A) > —2} if A e C e is 
&(\fe) close to \{ or 0(e) close to A + (0) or A~(0), so we only have to study A near these 3 points 
to determine the slow- fast eigenvalues of (|3.2| ). Note that the fast reduced (scalar) limit problem 
has an eigenvalue \! 2 = — | ((3.14), since Vo = 14(0) — > as e — * (Theorem 2.1)). We will prove 
below that £2 (A) has a (simple) zero close to A 2 , i.e. that the fast reduced eigenvalue A 2 persists. 

However, before going further into the details of the (possible) existence of eigenvalues near 
A 2 , A+(0) or A-(0), we formulate a result that is an immediate consequence of Lemma |4.2| and 
that establishes the stability of the wave for values of G\ and Hq such that the essential spectrum, 
and hence A + (0) and A~(0), is in the negative half-plane and not too close to the imaginary axis 
(see Lemma 3.2). 

Theorem 4.3. Let e > be small enough and let G\ < and H + G\ — 2t < be such 
that \G\\, \Hq + G\ — 2t| 3> e. The spectrum of the eigenvalue problem (3.1) associated to the 
stability of the solution e), V?,(£; e)) consists of a (unique) eigenvalue at A = and a part 

that is embedded in the region |Re (A) < — e}. Therefore, (Uh(£; s), Vh(£,; e)) is (s pect rally) stable. 

so that the 



0(v^ 



)■ 

We first note that indeed Af ast w^ = Q{eyfe) and A. s \ ovl v^ = 0(e), and 
for A E C e that are not asymptotically close to the possible degenerations 



Note that the operator defined by ( |3.l| ) is clearly sectorial in this case (see section 
nonlinear stability of (£/fo(£; e), T4(£; e)) follows by standard arguments (see for instance 

Proof of Lemma {A. We first note that indeed Af ast w^ = 0(e-\/e) and A s i 
thus < 2 (A) = 1 
of (|J) and (Q 

The inhomogeneous function u m may become unbounded as A approaches an eigenvalue, A{ = 
or A 2 = — |, or the essential spectrum al ss ( |3.15 ) of the linear problem associated to the fast 
reduced limit (3.12) with Vo = 0. To avoid irrelevant technicalities near cr^s we assume that 
A e C e n {Re(A) > — 2 + §}. The eigenfunction associated to A{, i.e. (£; 0) is odd, which implies 
that the inhomogeneous (and even) term Uh(£,;0) satisfies the solvability condition associated to 
(4.9) at A = 0. Hence, u ln remains bounded as A — > 0, so that £ 2 (A) = 1 + 0(^/e) also near A = 
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0, ||. The eigenfunction associated to X 2 is even, thus grows as 1/(A| — A) as A — > A2 H,|j|, 
which implies that i 2 (A,e) - 1 = 0{e^- a ) if |A - X f 2 \ = 0(e a ) for some a € (0 ^ 



' 2' 



The behavior of t2(X) near the degenerations of (4.7), i.e. the zeroes A (0) of Q(A; 0), follows 
from observing that A s i ow v^ = (t 2 — 1) x 0(e 1+<T ) if A is 0(e 2<T ) close to A + (0) or to A~(0) for some 

ae(o,|). ' □ 

It now follows, by a (standard) winding number argument jl], |], that the eigenvalue X 2 
persists as an eigenvalue of the full system ( J3 . l[ ) if it is not embedded in the essential spectrum. 

Lemma 4.4. Let G\ and Hq be such that a ess does not intersect an 0(e <T ) neighborhood of X 2 , 
for some a < ^. Then, there is an eigenvalue A2(e) of (3A) with lim £ ^o -^(e) = ^ 2 ~ ™ 



2 ' 



2 ■ 

Proof. By the assumptions in the Lemma, there exists a contour K in the complex A-plane 
that does not intersect cr css , that encircles an 0(e cr ) neighborhood of X 2 and that is 0(e CT ) close to 



X 2 . It follows from Lemma 4.2 that i2(A) = 1 + 0(e 2 <T ) for A e K , thus, the winding number of 



£2 (A) over if is 0. However, t 2 (X) must have a (simple) pole in the interior of K - as is observed 



in the proof of Lemma 4.2. We conclude that t 2 (X) must also have a (simple, real) eigenvalue in 
the interior of K. □ 
The possible existence of slow-fast eigenvalues near A + (0) or A~(0) is much more subtle. Since 
such eigenvalues only become relevant to the stability of the solution (Uh(£,; e), Vh(£; £ )) as G\ (or 



Ho + Gi — 2t) approaches (Theorem 4.3) we will consider this issue in the forthcoming sections. 



Remark 4.5. The eigenvalues \i(e) = and X 2 (e) — * — | as e — ► can be interpreted as 
'fast' eigenvalues, since they correspond to eigenvalues of the fast reduced limit problem. However, 
strictly speaking both eigenvalues also have the slow-fast structure described in the beginning of 
this section. 



First, we of course know that Ai(e) = is an eigenvalue - see also Lemma 3.8. Thus it is a 



zero of 2?(A, e). Since t 2 {X) = l + 0(y/e) for A near 0, see the proof of Lemma 4.2, we conclude that 
t\ (0; e) = (note that this in a sense obvious result does not follow directly from Lemma [T|). Thus, 
the solution 0i(£; 0, e) of (|3.2|) that by construction has a purely fast structure for £ <g; —1, does not 



blow up as e Al ^°' e ^ as £ — * 00 (Lemma 3.4). Nevertheless, the eigenfunction associated to A = 0, 



(^i,e(£)> ^h,i(0) nas a c l ear slow-fast structure, that it inherits from Vfc(£)) (Theorem 2.1). 

Hence, 0i(£; 0, e) is not the eigenfunction associated to A = 0. Neither is <fi 2 {£,', 0, e), since £2(0) ^ 0. 
It follows that the eigenfunction associated to A = must be a linear combination of (/>i(£;0,e) 
and (f> 2 (£,; 0, e), and thus that (/>i(£; 0, e) does not decay as £ — > 00, but instead grows linearly (and 
slowly), as e A2 ' 0,e ^ (like <p 2 (^;0,e)). The linear combination is such that the two growth terms 
e A 2 (o,e)C ^f or £ _> cancel. 

Second, X 2 (s) is not a zero of ti(A), although it is asympto tica lly close to such a zero, but it is a 



zero of *2(A). Thus, cj) 2 {^\ 2 {e),e) is the eigenfunction of (|3.2| ) at A = \ 2 {e) (and <px(£; A 2 (e), e) 
blows up fast, as e Al(A2(e) ' £)? ). 

4.2. The super-slow case: an example. In the previous section we have seen that the 
front might destabilize as G\ approaches (if we assume that Ho + G\ — 2r < 0). In this case, 
Theorem |2T] can no longer be used to establish the existence of the front (JJh{Q, Vh(Q)- Thus, 
the question about the stability of the front is closely related to the characteristics of the existence 
problem (as is usual in the analysis of (traveling) waves, see also Jl3||). In this section we consider 
the bifurcation as G\ approaches 0. Therefore, we assume that H — 2r < and 0(1) with respect 
to e. As in section || we consider in the super-slow case the simplified system in which the general 



function G(V) is replaced by a linear expression: G(V) — G{V — —e 2 ^V (see Remark 2.4). Note 
that Theorem |T3| a priori predicts a possible destabilization as G\ becomes 0(e), i.e. already 
before G\ — —7?, but it will be shown in the next section that the estimate in Theorem |4.3| is 
not sharp, in the sense that a bifurcation only occurs as G\ decreases to 0(e 2 ). 

One of the main differences between the analysis in this section and that of the regular case, 
is the fact that Vh(£) is no longer 0(e), i.e. V^(£) does not only contribute to the higher order 
terms in the stability analysis of the front solutions. Nevertheless, we follow the approach of the 
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previous section and express the solution u of (3.1) in terms of v, outside If (see (4.2)), 



Uu 



1 + V h - 3U 2 - A 



v + 0{e 2 u xx ) 



2(1 + Vh) + A 



■0(e 4 



v + 0(e 2 v xx ) 



(4.10) 



since 1 + T4(£;e) — C/^(£;e) = 0(e 4 ) (see fl2.8|) , recall that g 2 and G are 0(e 2 ) in the super-slow 
case). This yields that 

«« = {2 [i? (C/ h 2 , V},) + 0(e 4 )] C/fcW - [H (Ul V h ) + 0(e 4 ) - e 2 7 - rA] «} 



(4.11) 



= A 



H(l+V h ,V h ) 
2(l + V h )+\ 



+ e 2 1 + 0(e i )}v + 0(e 2 v xx ). 



It follows from section LI that one of the 'tips' of <t css , A + ( 0), i s 0(e 2 ) if G\ = 0(e 2 ) (and 
Ho — 2t < 0), while the other one, A~(0), is 0(1) and negative (|3.8|) . Thus, the destabilization of 
the front will either be ca used by a css at G\ = — 7, or possibly by a slow-fast eigenvalue A that 
is close to A + (0) (Lemma [O). Therefore, we introduce A by 



A = e 2 A, 



which implies that (4.11) can also be written as a super-slow system, 

H(i + v h ,v h y 



v xx = e < A 



2(1 + V h ) 



+ 7 + 0(e / ) }v 



(4.12) 



(4.13) 



As in section 2.2, we first consider the explicit example in which H(U 2 ,V) = HqU 2 . Thus, the 
existence of (several kinds of) front solutions is established by Theorem 2.3. In this case, the 
equation for v is, on the £-scale, given by 



A(r-iff o )+7 + 0(£ 2 ) 



v= [Al 3 (X,e) + 0(e 6 )]v, 



(4.14) 



see (3.4). Note that this equation is of constant coefficients type, and, at leading order, the same 
as in the equation for in the regular case (4.4). Hence, we can copy the arguments leading 
to Lemma LI a nd conclude that the fundamental solution $2 fe e 2 A, e) of (3.2) can again be 
expressed as in (4.5) outside the region If. Moreover, as in Lemma LI, we may conclude that 

i 2 (A,£)+i 3 (A,e) = l + 0(V^) ©■ 

We may now proceed as in the preceding section (and as in [Q, and determine £2 ( A) by 
measuring the change in the q = ^-coordinate of foiC) over the fast field. It follows from ( 4.14 ) 
that 



slow 



v 6 = 2e 2 (i 2 (A) - 1) A /A(r - -H ) + 7 + 0(e 2 V~e). 



(4.15) 



Note that we have to assume that A(r — |i?o) + 7 > 0, i.e. A 2 , 3 (A,e) > 0, which is a natural 
assumption, since 

2 7 



A 



tip 



A+(0) 



2t-H 



< 



(4.16) 



determines the 'tip' of a ess (recall that Ho — 2r < 0), i.e. £2 (A) is not defined if A < A t i p . By 
definition, Af aBt v^ is given by (4.8). Since, at leading order Vh(£,) = Vh{0) = v and Uh{£,) = 
u o(£,'t v o) (uniformly) in //(Theorem 2.3), and since uo((,;vo) decays exponentially fast on the 
(fast) £-scale, it follows that 

/>oo 

Afart^t = e 2 H / {2 [2u 2 {£; «o) - 1 - v ] u (£; v )u in (Z; v ) - ug(f; v )} d£ + 0(e 2 \/i), (4.17) 
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where Ui n (£;fo) is the (uniquely determined) bounded solution of the inhomogeneous problem 

Hi + (l + v - 3uq(£; v ))u = -Uo(f; v ). 

Since we already know one solution of the homogeneous problem, u(£) = uo,f (£; v o)j we can deter- 
mine u- m (^;vo) explicitly, 

Uin(C; v ) = ^77^ r (wo(C; vo) + £u(u(£; «o)) • (4.18) 

Thus, by ( |2.3| ), Af ast i>£ can be computed explicitly (at leading order), 

Afast^ = -e 2 H V2Vl + v + &(e 2 ^). 
Combining this with ( 4.15Q yields an explicit expression for t 2 {X), 



t 2 (X, e) = l- H 2(1 1 + J o) + 0(Vi) (4.19) 
V A(r-±ff ) + 7 



for A > A ti p ( |4.16| ). It follows that t 2 (X) > 1 + 0(y/e) for F < and t 2 {X) < 1 + 0(y£) for 



> 0. Hence, £2 (A) cannot have zeroes if Hq < 0. In other words, there cannot be an eigenvalue 



near the tip of the essential spectrum in case (ii) of Theorem 2.3. On the other hand, t 2 (X) can be 



for Hq > 0, i.e. in case (i) of Theorem 2.3 there indeed is a 'new' slow-fast eigenvalue of (3.1) 
it is given by 



Aedgc = e 2 A edgc = 27 + g o(l + ^o) £2 + 0{£2v - £) > = ^ (4 2Q) 

It — H 



( 4.16| ). Note that the eigenvalue A oc i gc merges with Ati p , and thus with <7 0SS as Hq J, 0. This 



is of course a leading order result, the accuracy of our analysis only allows us to conclude that 
I Atip — Aodgd < 0{e 2 y/e) as Hq J. 0, and that A e d g c does not exist for Hq < 0. Nevertheless, 
we conclude that A e dge appears from the essential spectrum as Hq increases through 0. In other 
words, A e dge is created, or annihilated, by an edge bifurcation. Note that the new eigenvalue 



appears exactly as <7 CSS becomes complex valued (section 3.1, Figure p.l| ). 

The existence, or non-existence of A e d g e is crucial to the character of the destabilization (see also 
the numerical simulations in section ||). For Hq < 0, the front solution (f/fc(£)> destabilizes 
as 7, or equivalently G\, crosses through 0. The destabilization is due to the essential spectrum, 
which implies that also the 'background states' (U(x,t),V(x,t)) = (±1,0) destabilize at 7 = 0. 



However, in the case Hq > the eigenvalue is A e d ge is e 2 Hq{\ + v )/ (2t — Hq) ahead of cr CS s ( 4.20| ) 



in the sense that it reaches the axis Re(A) = before fj C ss as 7 > decreases to 0. Thus, if Hq > 
the front solution (Uh(^),Vh{£,)) destabilizes by an element of the discrete spectrum of (|3.l|) at 
7doubio = \Hq{1 + v ) + 0(v / £ : ) > 0. As a consequence, the background states (±1,0) remain 
stable as (Uh(£), VJi(£)) destabilizes for Hq > 0, contrary to the case Hq < 0. The bifurcation at 



7doubie is associated to the saddle-node bifurcation of heteroclinic orbits described in Theorem 2.3 



Theorem 4.6. Assume that G(V) = -s 2 j, H(U 2 ,V) = H U 2 , Hq - 2r < and 0(1), and 
that e > is small enough. 

1 ^ 1 _i_ ^ _i_ 2 

(i) Let (U h ' (£),V h ' (£)) and (U h ' (£),V h ' (£)) be the two types of heteroclinic front solutions 
that exist for Hq > and 7 > 7d ou bio = S-ETq + 0(V^) with, at leading order, < V^~' (0) = v% < 
2 < v 2 = V /l +,2 (0) (Theorem \2.!\ ). The front solution (U^' 1 ^), Vf^' 1 ^)) is asymptotically stable for 

1 ^ 1 2 _i_ ~^ 1 

7 > 7doubio> the front (U h ' (£),V h ' (£)) unstable; (Uj[' (£),V h ' (£)) destabilizes by an element of 
the discrete spectrum, A c dg C ; at 7 = 7doubic o,nd merges with (U^' (£), (£)) in a saddle-node 
bifurcation of heteroclinic orbits. 
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(ii) Let (Ujt(£), Vh~(tj)) be a heteroclinic front solution that exist for Hq < and (all) 7 > 



(Theorem 2.t ); (C/^~(£), V^ + (£)) is asymptotically stable for all 7 > 0, it is destabilized at 7 = by 
the essential spectrum o- css . 

Remark 4.7. As in the regular case, spectral stability implies asymptotic nonlinear stability in 
this super-slow case, since the linear operator associated to the stability problem remains sectorial 
as long as e > 0. 

Proof of Theorem (.6. We first note that the condition Ho — 2t < an d (1) determines that 
er css can only cross, or come close to, the Re(A) = 0-axis at A = (Lemma 3.2 with G\ = 0(e 2 )). 



(i) The eigenvalue 'in front of er css , \^ s<} {vi,2), is given by (4.20), where vq > is a solution of 
97V 2 = 2H 2 (1 + v) 3 , and v = v 1 < 2 (at leading order) for (U^' 1 ^), while v = v 2 > 2 
(at leading order) for (U£' 2 (£), V^' 2 ^)) - Theorem |J. Thus, by ( pL2C| ) — A^ dgc (^i) < and 
A 2 dgc (w2) > if 7 < 7doubio = §#0 + 0(\/e)- As a consequence, Aj dgc (ui) | and A 2 dgc (u 2 ) I as 
7 I 7doubic, at which the saddle-node bifurcation takes place.. 

(ii) We have already shown that there can be no eigenvalues in front of the tip of cr ess . Therefore, 
the statement of the Theorem follows. □ 

Remark 4.8. Since i 2 (A) = 0, the slow-fast eige nfun ction associated to the bifurcation at 
7 = 7doubie is given by 4>2(£,)- It follows from Lemmas 3/7 and 3^? that the u-component of </>2 is 
odd, and the u-component even, as function of £. 

4.3. Bifurcations in the general super-slow problem. We now consider the stability of 
a front solution in the general super-slow limit. Thus, we assume we have est ablis hed the existence 
of a front (Uh(£,), Vh{C)) f° r a certain given function H(U 2 ,V) (Theorem 2J5). To analyze its 
stability, we again try to determine £2 (A) by measuring Af as tU|_and A s i ow ^. 

In order to determine A s i ow t>£ we follow the derivation of ( 4.13 ) in the previous section. Hence, 
we again con clude that non-trivial eigenvalues near are only possible for A = 0(e 2 ), thus we again 
introduce A ( 4.12| ) (see also the proof of Theorem 4.1C for more details on the necessity of this 
scaling). Note that both G\ and A are now 0(e 2 ), thus, we can immediately obtain a leading order 
expression for Af ast w^ in terms of H(U 2 , V), 



Af ast l>£ 



; n o {2[i?(u 2 , W0 )-(i+ 



Vq 



u o)§m( u o, v o)] u u h 



- [H(u 2 Q , vo) - (1 + «o - u 2 )^ (u 2 , v Q )] } d£ + 0(e 2 V^) 



(4.21) 



(3.1), where Mi n (£) is given in ( 4.18 ) - rec all t hat v — 1 + &(^fe) in //. As in the p revi ous section, 
we have approximated Uh{£) by u Q (^; vq) (|2.3|), Vh(£) by v and If by M (Theorem 2.5). Note that 
the integral converges and that Af ast Vf is (at leading order) independent of 7 and A. 

It is in principle possible to determine A s i ow w^ in terms of £2 (A) from (4.13), however, this 
equation is in general not of constant coefficients type (unlike for the example problem in section 
p~2] ). If we introduce the super-slow coordinate X by X = ex = e 2 £, we can write ( 1.13 ) as 



vxx 



A 



H(l + V h (X),V h (X)) 
2{l + V h {X)) 



+ 7 + 0(£ 2 ) 



(4.22) 



i.e. the functions Vh(X) introduce explicit X-dependent terms in the equation (section |2T3| , Vh(X) 
behaves as e T ^ x on -Mf). Nevertheless, we can in principle determine the u-components of the 
solution 02(C) °f (3.2) outside th e fas t region If. H owev er, the analysis is much less transparent. 
For instance, the decomposition ( |4.5| ) as in Lemma 11 now only holds for X ^> 1, therefore the 
relation between £3 (A) and £2 (A) that is obtained from the value of v in If will in general be more 

H(l + V h (X),V h (X)) ~ 
2(l + V h (X)) 



complicated than in (4.6). Moreover, A 



7 might change sign as function 

of X, so that the solution v(X) of Q4.22) ) can have oscillatory parts. 

Thus, we conclude that it is not a straightforward extension of the approach in previous section 
to determine *2(A) for general values of A. It should also be noted that a similar problem occurs in 
the regular case, in the study of possible eigenvalues near A ± (0) (Lemma 4.2). If one introduces 
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A ± by A = A ± (0) + eA*, a nd de rives the leading order equation for v xx (4.2) in this case, then 
one finds an equation like (4.22), i.e. an equation with spatially dependent coefficients (these 
rr-dependent terms originate from the 0(e) corrections corresponding to Vh(x) = 0(e) in (4.2) 
and ( |4.3| )). Hence, at this point it is not yet possible to determine in full detail whether or not 
eigenvalues exist near the tips of a ess for general nonlinearities H(U 2 , V) and general A. Moreover, 
it is also not possible to explicitly describe how and when eigenvalues appear from, or disappear 
into, (Joss- On the other hand, it is clear from (4.21) and (4.22) that the number of zeros of ^(A) 
depends (for instance) on Hq. It thus follows that eigenvalues will be created/annihilated near the 
tip of (Jess in the general case (like in the example system considered in the previous section) . The 
analysis of eigenvalues near the tip of o~ ess is therefore a continuing subject of research (in progress, 
see also section ||) . 

Nevertheless, the val ue A = A = is, of course, especially relevant for the stability analysis of 
the front, and, equation ( [4.22] ) is again of constant coefficients type at lea ding order for this special 
value of A. Hence, for A = we can obtain the equivalent of Lemma 4.1, so that it follows that 



A slow z; ? U = o = 2e 2 (i 2 (0) - 1)^7 + 0(e\/i). 



(4.23) 



Note that eventually it becomes clear at this point why the choice G\ = — e 2 7 is the most relevant 
scaling of G±. With this scaling the 'jumps' A s i ow w^ and Af ast i>£ (4.21) are of the same magnitude 
in e at A = 0. Therefore, £2(0, e) is asymptotically close to 1 for all G\ with |Gi| 3> e 2 - see 
Lemma [4.2| and its proof. Thus, the stability problem (3.1) can only have a double eigenvalue at 
if G\ = 0(e 2 ). This establishes a significant link between the stability analysis and the existence 
analysis of section g, since it is clear from the analysis there that the sc alin g G± = 0(e 2 ) is also 
the most relevant scaling for the (super-slow) existence problem (Remark Moreover, this link 
is even much more explicit. 

Theorem 4.9. Assume that G(V) = — e 2 j, that H — 2t < and 0(1), and that e > is 
small enough. Let the front solution (Uh(£,; e), Vjj(£; e)) be a heteroclinic solution that corresponds 
to an intersection T~ H W u (— 1, 0, 0, 0)1^- as described in Theorem 2.1. The stability problem 
associated to the front solution has a double eigenvalue at X — if and only if the intersection 
T~ C\W U (— 1, 0, 0, 0)1^- is non-transversal. If the intersection T~ C\W U (— 1,0,0, 0)\ M - is a second 
order contact, then the front bifurcates at 



< 7doublc — 



_ 1 



4(l + u ) 



Cjl + va-uDHiulvaW 



(1 



«o)f£««o)R 



(4.24) 



by merging with another front solution in a saddle-node bifurcation of heteroclinic orbits. 

Proof. First, we recall from section 2.3 that a heteroclinic connection tha t cor responds to 
the intersection of ^"(—1,0,0,0)1^- = {q = £t/jv} and T~ is determined by ( |2.15| ). This is of 
course a leading order approximation. In the proof of this Theorem we refrain from mentioning 
this obvious fact at several places. To determine the ^o-dependence of the right hand side of this 
relation, we define u>o(£) as the (monotonically increasing) heteroclinic solution of w+(\— w 2 )w = 0. 
It follows that 



uo(£; w o) = Vl + wu>o(vl + vo£), w (t) = t&nh\/ -t 



(4.25) 



(2.3). Replacing u (^;v Q ) by w Q (t) in ( [2.15] ) yields 

Vivo = -v 7 ! + vo J (l-w )H((l + v )w Q ,v )dt. 



(4.26) 
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Thus, T~ n ^"(-1,0,0,0)1^- is non-transversal if (flf) holds and 

\fr = hmr {VTT^J^a - w 2 )H((i + v )w 2 ,v )dt} 
= j 7 ^C oo (i-<)H((i + v )wlv )dt 



+ iyiT^/_ 00 00 (l - wl)[wl^{(\ + v )w 2 ,v Q ) + ((1 + v )w 2 ,v )]dt (4-27) 



OH , 

+ «o - «§)[«§|^(«§, «o) + (i + «d)f£(«§, voM, 



by re-introducing mq(£; up ). No te th at (4.24 ) follows from this equation. The expression for £2(0, 
is determined by ( |4.2l| ), ( jl.23| ) and (|4.18|) , 



t 2 (0,e) = l- 



2V7(l + «o) 



0(v^), 



where 



Zi 

1 3 



J-oot 1 + v o- t4)H(v%, v )d£, 
o)l a o dl 

f-ocK 1 + v ° ~ u a)m?( u i v o) ~ H(ulv )]^u u ^. 



= .TJI + ^0 - u 2 )[u 2 §^(u 2 ,v ) + (1 + Wo )ff ( Uo > )K, 



(4.28) 



We find by partial integration that 
Z 3 = 



-Ji, 



which implies that 



t 2 (0,e) = 1- 



2Z, 



4^7(1 + ^0) 



so tha t we can conclude by ( 4.28 ) that ^(0, e) = is equivalent to the non-transversality condition 
(4.27). Hence, a double eigenvalue of (3.1) coincides with a saddle-node bifurcation of heteroclinic 
orbits, unless the tangency between T~ and ^"(—1,0,0,0)1^- is degenerate. □ 
Finally, we can turn to the question about the char acte r of the destabilization of the regular front 
solution, that has been studied in sections 2.1 and 4.1, as G\ approaches 0. In order to do so, we 
first note that the existence problem for the regular case can be recovered from that of the singular 
limit by re-introducing G\ = — 'ye 2 in the existence condition (2.15). This implies that vo must 
become 0(e) and that 



1 r°° 

\f Z Glv Q = e- (1 - ul)H(ul OR 
* J — 00 



(4.29) 



which is equivalent to (2.4) in Theorem pT[ Thus, the structure of the front (Uh(£), Vh(£)) as 
function of G\ f can be determined by tracing the intersection T~ n W u (— 1, 0, 0, 0)|^- in 
the super-slow limit as W u (— 1,0,0, 0)\ M - — {q = e^jv} goes down from being almost vertical 
(Gi = 0(1), 7 = 0(l/e 2 )) to horizontal (G\ =7 = 0). Note that this process determines a 
unique 'regular' element in the intersection T~ n W u {— 1,0,0, 0)\ M - , all other elements of T~ n 
W"(— 1, 0,0, 0)1^- do not persist in the regular limit 7 = 0(l/e 2 ) (here, we do not pay attention 
to possible heteroclinic connections that have Do > 1 as 7 > 1). I t depends on the sign of 
\ (1 ~ u o)H(u 2 ) ,0)dt; whether vq will be positive or negative ( 4.29| ), i.e. whether the regular 
intersection T~ n W u (— 1,0,0, 0)\ M - travels through the first or through the third quadrant of 
the (v,<7)-plane as 7 decreases. Since H(U 2 , V) is smooth, we can make a distinction between two 
different types of behavior: 
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Type D The regular element of T~ n W u {— 1, 0, 0, 0)1^- merges at a certain critical value of G\ = 
—e 2 j < with another element of T~ n W u (— 1, 0, 0, 0)| w- in a saddle-node bifurcation 
of heteroclinic orbits. 

Type E The regular element of T~ n W u (-1, 0, 0, 0)1^- exists up to the limit G\ = 0. 
Note that T~ approaches (-1,0) as v J, -1 ( |4.26| ), so that an element of T~ n 0, 0, 0)\ M - 

can only reach the singular region {vq < —1} at 7 = 0, which indeed implies that there can only 
be orbits of type D and E in the third quadrant. We can now describe the destabilization of the 
regular fronts as G\ approaches 0. 

Theorem 4.10. Assume that G(V) = -£ 2 7, that H - 2r < and 0(1), and that e > is 
small enough. Consider the heteroclinic front solution (f/ft(£), Vh(£)) determined in Theorem 2.1 



for G\ < and 0(1) and in Theorem \2. j| for G\ — 0(e 2 ). If the fr ont is of type D as G\ becomes 
0(e 2 ), then it is asymptotically stable up to G\ — — £ 2 7doubic < ( 4-^4 ) an d it is destabilized by 
a (discrete) eigenvalue through a saddle-node bifurcation of heteroclinic orbits. A front solution of 
type E is stable up to G\ = and it is destabilized by the essential spectrum. 

Thus, the destabilization of a regular front solution in the limit G\ | is completely determined 
by the geometrical structure of T~ n W u {— 1, 0, 0, 0)1^- in the super-slow limit. Note that Figure 
[2~5| presents examples of type D and type E behavior. 

Proof. The proof of this Theorem is a bit more subtle than a priori might be expected, 
since in general we do not have control over the eigenvalues of (3.1) near the tip of er ess (see also 
Remark 4.11), except that these eigenvalues must be 0(e 2 ) close to a css (see also below). Thus, 
for instance the following scenario for a type D orbit might be possible as 7 decreases to 7doubi e : 
two eigenvalues bifurcate (subsequently) from cr ess (as real eigenvalues), merge and become a pair 
of complex eigenvalues. This pair crosses through the Re(A) = axis at 7Hopf > 7dout>ie, and touch 
down again on the real axis. At 7doubic one of these eigenvalues returns to Re(A) = 0. Thus, 
in this scenario, there already exists an unstable eigenvalue at 7 — 7doubio; moreover, the front 
destabilizes by a Hopf bifurcation at 7Hopf > 7doubie- 

Let us first note that a destabilization by a Hopf bifurcation is the only alternative to the 
statements of the Theorem, since eigenvalues either move through 0, or (in pairs) through the 
Re(A) = axis. If we can show that a Hopf bifurcation cannot occur for 7 > 7<j ubie; then it is 
clear that for type D orbits A c d gc < for 7 > 7doubie and that there is no unstable spectrum at 



7 = 7doubie (this follows from Theorem 4.3: if 7 is ^> 0(l/e), all non-trivial eigenvalues must be 
in {Rc(A) < — e}, hence, by decreasing 7, there is one e igen value, A d ge , that is the first to reach 0; 
this happens at the saddle-node bifurcation (Theorem |4.9| ), i.e. at 7 = 7doubic)- Thus, the front is 
stable for 7 > 7doubic- The same argument can be used to establish the non-existence of unstable 
spectrum for type E orbits, if there are no Hopf bifurcations possible. 

To show that there cannot be Hopf bifurc ation s (for Ho — 2r < and 0(1), see section we 
first ascertain that A must be 0(e 2 ), i.e. that ( 1.12| ) is the correct scaling. This follows by the same 
arguments as in the proof of Lemma 4.2. If A s i ow i>£ 3> Af ast i>£, then there cannot be an eigenvalue. 
Thus, it follows from (4.11) that |A| must indeed be 0(e 2 ) near A + (0). Hence, even if there is 
a Hopf bifurcation, it will be 0(e 2 ) close to 0. Next, we realize that this situation is covered by 
( 4.17 ) for the jump through the fast field, thus, Af ast t>£ is real (at leading order), independent of 
A. However, it follows from (4.22) that A s i ow U£ cannot be real if A is complex- valued. Hence, there 
cannot be a Hopf bifurcation 0(e 2 ) close to A = 0. □ 

Remark 4.11. By the same geometrical arguments (that are based on Theorem 4.9) we can 
describe the character of the bifurcations as function 7 in the stability problem associated to a 
heteroclinic orbit that corresponds to a non-regular element of T~ n W u (— 1, 0, 0, 0)1^- . However, 
it should be noted that, in general, we do not have enough information on the spectrum of ( |3.l| ) 
to establish the stability of such a front, since we did not determine all possible eigenvalues. 
In general, we cannot exclude the possibility that various eigenvalues have bifurcated from the 
essential spectrum for these fronts (in fact, the possible oscillatory character of a solution v(X) of 
( 4.22j ) strongly suggests that this can happen). Nevertheless, we may for instance conclude that if 
the regular orbit is of type D, then it merges with a non-regular orbit at 7doubie that is unstable 
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for any 7 > 7doubie for which it exists. 

Remark 4.12. The most simple example one c an c onsider is H(U 2 ,V) = Ho- This corre- 
sponds to the case in which the function F(U 2 , V) in (1.1) is (the most general) linear function of 
U 2 and V with parameters d and H (i.e. F(U 2 ,V) = F + (# +Gi)V - H U 2 , recall that F(l, 0) 
must be 0). In this case, T~ is given by {q = 2eH ^l + v + 0(e 2 )}, so that W U (-1,0,0,0)\ M - 
can never be tangent to T~ . Hence in this case, there is a uniquely determined front solution 
of type E for any Ho ^ and G\ < 0, i.e. the front solution is stable up to G\ = and it is 
destabilized by the essential spectrum. 

Remark 4.13. We did not consider the degenerate case in which H(U 2 ,V) is such that 
H(l + V, V) = (section @), i.e. functions H such that H{U 2 , V) = (1 + V - U 2 )H(U 2 , V) for 
some smooth function H. In a sense, this is a much more simple problem, for instance since in 
the super-slow limit, the stability problem in the slow field is automatically of constant coefficients 



type (at leading order), see (4.11), (4.22). Moreover, it is also clear from these same relations that 
we can find 0(1) instead of 0(e 2 ) eigenvalues in this case if r = 0(e 2 ). In fact, the situation 
is very much like the stability analysis of (homoclinic) pulses in mono-stable systems in j|, 
For instance, as in |5|, potential eigenvalues are no longer 'slaved' to the tips of the essential 
spectrum or to the eigenvalues of the fast reduced limit (Lemma 4.2). Moreover, the 'natural' 
persistence result of Lemma 4.4 is also not valid in this case, in general. 



5. Simulations and Discussion. 

5.1. Simulations. We now e xam ine nume rically the difference between the two types of 
bifu rcat ions discussed in Theorems 4.6 and 4.1C . We consider the example system of sections 2.2 
and 4.2 for Ho > (case (i), type D) and Ho < (case (iiO, type E). First, we note that in 
both cases the simulations confirm that the fronts are asymptotically stable up to the analytically 
determined bifurcation values. In case (i) the front destabilizes at 7 < 7dout>ie due to an eigenvalue 
in the discrete spectrum. The eigenfunction associated to this type of destabilization is localized 
to a neighborhood of the front as can be seen in figure 5.1. In this case the front becomes unstable 
and blows up in finite time, while the background states remain stable. In case (ii), the tip of 
the essential spectrum becomes positive and the background states become unstable as 7 passes 
through 0. As can be seen in figure 5.1, this destabilization causes the front to collapse. The 
U component then tends to on the entire real line and the V component grows according to 
Vt — V xx +£ 2 \"(\V . Thus, we may conclude that type D or type E orbits indeed exhibit significantly 
different behavior at the destabilization. These simulations were performed using SPMDF 
with Neumann boundary conditions at x = ±50. The initial conditions used in figure (5.1) are 

/Mx| 



given by U (x, 0) = uq(x, e; vi) (2.3 and V(x, 0) = V\e 



(as described in Theorem 2.3) 



5.2. Hopf bifurcations. As we have seen in section 4.1, in general there can be (complex) 
eigenvalues near the endpoints A ± (0) of a CBS . Thus, if we keep G\ < fixed at an 0(1) value and 
increase Ho such that Ho + G\ — 2t approaches 0, we encounter a similar issue as was studied in 
the previous section: will the front be destabilized by a css at Ho — It — G\ or (just) before that, 
by an eigenvalue? In this case, the bifurcation is of Hopf type, and it is not associated to the 
existence problem. This problem can in principle be analyzed by the methods developed here, i.e. 
by determining t 2 (X, s) through A s i ow U£ and Af ast f£-_We have already mentioned the new features 
of the measuring the slow 'jump' A slow v 5 in section L3. Moreover, since the bifurcation does not 
occur near A = 0, we do not have an explicit formula for Uin(£)i like (4.18), and it is thus not 
immediately clear whether it is possible to determine Af ast V£. Note that this latter issue is solvable 
with the hypergeometric functions method developed in j| ||. Nevertheless, we do not go deeper 
into this subject here. 

5.3. Planar fronts and stripes. A next step in the study of (planar) stripes, as mentioned 
in the Introduction, is the stability analysis of planar fronts, i.e. the analysis of the stability 
of the fronts (Uh(Q,Vh(£)) with respect to two-dimensional perturbations (thus, (Uh(£),Vh(£)) 
represents a planar front that has a trivial structure in the y-direction) . The methods developed 
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here can be used to study this problem (as is also suggested by in which a similar problem has 
been studied in a mono-stable Gierer-Meinhardt context). It should be noted here that there are 
several papers in the literature that consider the question of the (non-)persistence of the stability 
of one-dimensional fronts as two-dimensional planar fronts (see for instance |l(| [H^, [l^, |l5|). The 
analysis in [jl9| [l^] of a class of singularly perturbed bi-stable systems shows that the planar fronts 
considered there cannot be stable, while it is shown that planar fronts can be stable in a more 
regular context in Jl2[ ] . Thus, this is a nontrivial issue. Preliminary analysis of the front solutions 
considered in this paper indicates that these solutions remain stable as planar fronts in the regular 
case (i.e. as long as G\ < and 0(1)). The analysis of the planar fronts, and their spatially 
periodic counterparts, the stripe patterns, is the subject of work in progress. 
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doctoral fellowship. A.D. and D.I acknowledge support of the 'Research Training Network (RTN): 
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(a) Re(A) vs k with H < 0. 



(b) Re(A) vs Im(A) with H < 0. 




(i) H € ({V^f- v /r GT) 2 ,2r- Gi). 



(j) ffo S ((v 7 ^- v /r GT) 2 ,2r-Gi). 



Fig. 3.1. The five possible different structures of the stable essential spectrum. On the left we plot Re(X) vs 
on the right Re{\) vs Im{\). 
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(a) (/-component 



(b) V-componcnt 



Fig. 5.1. Numerical simulation of destabilization caused by the discrete spectrum, both components blow up in 
finite time (Ho = — 1, 7 = 1.4, r = 1 and e = O.lj. 



-0.2 - 
-0.4 - 
-0.6 - 



(a) {/-component 



(b) V-componcnt 



Fig. 5.2. Numerical simulation of destabilization caused by the essential spectrum, U — > and \V\ grows 
slowly and exponentially (Ho = I, 7 = —0.1, r = 1 and e = O.lj. 
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Abstract. Documentation is given for use of the SIAM IATf^X macros. These macros are now 
compatible with I^T£]X2 e . Instructions and suggestions for compliance with SIAM style standards 
are also included. Familiarity with standard IATjrjX commands is assumed. 

Key words. 

AMS subject classifications. 

1. Introduction. This file is documentation for the SIAM DTjtX macros, and 
provides instruction for submission of your files. 

To accommodate authors who electronically typeset their manuscripts, SIAM sup- 
ports the use of IATfrX. To ensure quality typesetting according to SIAM style stan- 
dards, SIAM provides a I^Tj^X macro style file. Using to format a manuscript 
should simplify the editorial process and lessen the author's proofreading burden. 
However, it is still necessary to proofread the galley proofs with care. 

Electronic files should not be submitted until the paper has been accepted, and 
then not until requested to do so by someone in the SIAM office. Once an article is 
slated for an issue, someone from the SIAM office will contact the author about any 
or all of the following: editorial and stylistic queries, supplying the source files (and 
any supplementary macros) for the properly formatted article, and handling figures. 

When submitting electronic files (electronic submissions) (to tex@siam.org) in- 
clude the journal, issue, and author's name in the subject line of the message. Authors 
are responsible for ensuring that the paper generated from the source files exactly 
matches the paper that was accepted for publication by the review editor. If it does 
not, information on how it differs should be indicated in the transmission of the file. 
When submitting a file, please be sure to include any additional macros (other than 
those provided by SIAM) that will be needed to run the paper. 

SIAM uses MS-DOS-based computers for IATj^X processing. Therefore all file- 
names should be restricted to eight characters or less, plus a three character extension. 

Once the files are corrected here at SIAM, we will mail the revised proofs to be 
read against the original edited hardcopy manuscript. We are not set up to shuttle 
back and forth varying electronic versions of each paper, so we must rely on hard copy 
of the galleys. The author's proofreading is an important but easily overlooked step. 
Even if SIAM were not to introduce a single editorial change into your manuscript, 
there would still be a need to check, because electronic transmission can introduce 
errors. 

The distribution contains the following items: siamltex . els, the main macro 
package based on article. els; siamlO . clo, for the ten-point size option; 
subeqn. clo, a style option for equation numbering (see §3 for an explanation); and 
siam.bst, the style file for use with BibTjtX. Also included are this file docultex.tex 
and a sample file lexample.tex. The sample file represents a standard application 
of the macros. The rest of this paper will highlight some keys to effective macro use, 
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as well as point out options and special cases, and describe SIAM style standards to 
which authors should conform. 

2. Headings. The top matter of a journal paper falls into a standard format. 
It begins of course with the \documentclass command 

\documentclass{siamltex} 

Other class options can be included in the bracketed argument of the command, 

separated by commas. 

Optional arguments include: 

final Without this option, lines which extend past the margin will have black boxes 
next to them to help authors identify lines that they need to fix, by re- writing 
or inserting breaks, final turns these boxes off, so that very small margin 
breaks which are not noticible will not cause boxes to be generated. 

oneeqnum Normally siamltex. els numbers equations, tables, figures, and theorem 
environments with a decimal number, composed of the section of the paper, 
a period, and the number of the enumerated object (example: 1.2). The 
sequence of numbering is also restarted with each new section, so that, for 
example, the last equation of section 3 may be 3.10, but the first equation of 
section 4 would be 4.1. Using oneeqnum numbers all equations consecutively 
throughout a paper with a single digit. 

onethmnum Using onethmnum numbers all theorem-like environments consecutively 
throughout a paper with a single digit. 

onefignum Using onethmnum numbers all figures consecutively throughout a paper 
with a single digit. 

onetabnum Using onethmnum numbers all tables consecutively throughout a paper 
with a single digit. 

The title and author parts are formatted using the \title and \author com- 
mands as described in Lamport ||. The \date command is not used, \maketitle 
produces the actual output of the commands. 

The addresses and support acknowledgments are put into the \author commands 
via \thanks. If support is overall for the authors, the support acknowledgment should 
be put in a \thanks command in the \title. Specific support should go following 
the addresses of the individual authors in the same \thanks command. 

Sometimes authors have support or addresses in common which necessitates hav- 
ing multiple \thanks commands for each author. Unfortunately I^Tj^X does not nor- 
mally allow this, so a special procedure must be used. An example of this procedure 
follows. Grant information can also be run into both authors' footnotes. 

\title{TITLE OF PAPER} 

\author{A. ~U. Thorone\f ootnotemark [2] \ \f ootnotemark [5] 
\and A.~U. Thortwo\f ootnotemark [3] \ \f ootnotemark [5] 
\and A.~U. Thorthree\f ootnotemark [4] } 

\begin{document} 
\maketitle 

\renewcommand{\thef ootnote}{\f nsymbol{f ootnote}} 
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\f ootnotetext [2] {Address of A.~U. Thorone} 

\footnotetext [3] {Address of A.~U. Thortwo} 

\f ootnotetext [4] {Address of A.~U. Thorthree} 

\f ootnotetext [5] {Support in common for the first and second 

authors . } 

\renewcommand{\thef ootnote}{\arabic{f ootnote}} 

Notice that the footnote marks begin with [2] because the first mark (the as- 
terisk) will be used in the title for date-received information by SIAM, even if not 
already used for support data. This is just one example; other situations follow a 
similar pattern. 

Following the author and title is the abstract, key words listing, and AMS subject 
classification number(s), designated using the {abstract}, {keywords}, and {AMS} 
environments. If there is only one AMS number, the commands \begin{AM} and 
\end{AM} are used instead of {AMS}. This causes the heading to be in the singular. 
Authors are responsible for providing AMS numbers. They can be found in the 
Annual Index of Math Reviews, or through e-Math (telnet e-math . ams . com; login 
and password are both e-math). 

Left and right running heads should be provided in the following way 

\pagestyle{myheadings} 
\thispagestyle{plain} 

\markboth{A.~U. THORONE AND A.~U. THORTWOMSHORTER PAPER 
TITLE} 

3. Equations and mathematics. One advantage of P/TjtX is that it can auto- 
matically number equations and refer to these equation numbers in text. While plain 
Tj^X's method of equation numbering (explicit numbering using \leqno) works in the 
SIAM macro, it is not preferred except in certain cases. SIAM style guidelines call for 
aligned equations in many circumstances, and P>TjtX's {eqnarray} environment is not 
compatible with \leqno and P>TjtX is not compatible with the plain TjtX command 
\eqalign and \leqalignno. Since SIAM may have to alter or realign certain groups 
of equations, it is necessary to use the I^TjtX system of automatic numbering. 

Sometimes it is desirable to designate subequations of a larger equation number. 
The subequations are designated with (roman font) letters appended after the number. 
SIAM has supplemented its macros with the subeqn.clo option which defines the 
environment {subequations}. 

\begin{subequations}\label{EKx} 
\begin{equation} 

y_k = B y_{k-l} + f , \qquad k=l,2,3,\ldots 
\end{equation} 

for any initial vector $ y_0$. Then 
\begin{equation} 

y_k\rightarrow u \mbox{\quad iff\quad} \rho( B)<1. 
\end{equation} 
\end{subequations} 
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All equations within the {subequations} environment will keep the same overall 
number, but the letter designation will increase. 

Clear equation formatting using T}rX can be challenging. Aside from the regular 
TgX documentation, authors will find Nicholas J. Higham's book Handbook of Writing 
for the Mathematical Sciences useful for guidelines and tips on formatting with 
TgX. The book covers many other topics related to article writing as well. 

Authors commonly make mistakes by using <, >, \mid, and \parallel as delim- 
iters, instead of \langle, \rangle, I , and \ I . The incorrect symbols have particular 
meanings distinct from the correct ones and should not be confused. 



Table 3.1 
Illustration of incorrect delimiter use. 



Wrong 


Right 


<x, y> < x, y > 
5 < \mid A \mid 5 <| A \ 
6x = \parallel x 

- l\parallel_{i} Qx =| x — 1 \\i 


Mangle x, y\rangle (x,y) 
5 < |A| 5 < \A\ 

6x = \|x - l\|_{i} 6x = \\x - l||j 



Another common author error is to put large (and even medium sized) matrices 
in-line with the text, rather than displaying them. This creates unattractive line 
spacing problems, and should be assiduously avoided. Text-sized matrices (like (^ b )) 
might be used but anything much more complex than the example cited will not be 
easy to read and should be displayed. 

More information on the formatting of equations and aligned equations is found 
in Lamport ||. Authors bear primary responsibility for formatting their equations 
within margins and in an aesthetically pleasing and informative manner. 

The SIAM macros include additional roman math words, or "log- like" functions, 
to those provided in standard TgX. The following commands are added: \const, 
\diag, \grad, \Range, \rank, and \supp. These commands produce the same word 
as the command name in math mode, in upright type. 

4. Special fonts. SIAM supports the use of the AMS-TgX fonts (version 2.0 and 
later). The package amsfonts can be included with the command 
\usepackage{amsf onts}. This package is part of the AMS-I^TgXdistribution, avail- 
able from the AMS or from the Comprehensive TeX Archive Network (anonymous 
ftp to ftp.shsu.edu). The blackboard bold font in this font package can be used for 
designating number sets. This is preferable to other methods of combining letters 
(such as I and R for the real numbers) to produce pseudo-bold letters but this is 
tolerable as well. Typographically speaking, number sets may simply be designated 
using regular bold letters; the blackboard bold typeface was designed to fulfil a desire 
to simulate the limitations of a chalk board in printed type. 

4.1. Punctuation. All standard punctuation and all numerals should be set in 
roman type (upright) even within italic text. The only exceptions are periods and 
commas. They may be set to match the surrounding text. 

References to sections should use the symbol §, generated by \S. (If the reference 
begins a sentence, the term "Section" should be spelled out in full.) Authors should 
not redefine \S, say, to be a calligraphic S, because \S must be reserved for use as the 
section symbol. 
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Authors sometimes confuse the use of various types of dashes. Hyphens (-, -) are 
used for some compound words (many such words should have no hyphen but must 
be run together, like "nonzero," or split apart, like "well defined"). Minus signs ($-$, 
— ) should be used in math to represent subtraction or negative numbers. En dashes 
( — , -) are used for ranges (like 3-5, June- August) , or for joined names (like Runge- 

Kutta). Em dashes ( , — ) are used to set off a clause — such as this one — from the 

rest of the sentence. 

4.2. Text formatting. SIAM style preferences do not make regular use of the 
{enumerate} and {itemize} environments. Instead, siamltex.cls includes defini- 
tions of two alternate list environments, {remunerate} and {romannum}. Unlike the 
standard itemized lists, these environments do not indent the secondary lines of text. 
The labels, whether defaults or the optional user-defined, are always aligned flush 
right. 

The {remunerate} environment consecutively numbers each item with an arabic 
numeral followed by a period. This number is always upright, even in slanted envi- 
ronments. (For those wondering at the unusual naming of this environment, it comes 
from Seroul and Levy's 0] definition of a similar macro for plain TpjX: \meti which 
is \item spelled backwards. Thus {remunerate} a portion of {enumerate} spelled 
backwards.) 

The {romannum} environment consecutively numbers each item with a lower-case 
roman numeral enclosed in parentheses. This number will always be upright within 
slanted environments (as in theorems). 

5. Theorems and Lemmas. Theorems, lemmas, corollaries, definitions, and 
propositions are covered in the SIAM macros by the theorem-environments {theorem}, 
{lemma}, {corollary}, {definition} and {proposition}. These are all numbered 
in the same sequence and produce labels in small caps with an italic body. Other 
environments may be specified by the \newtheorem command. SIAM's style is for 
Remarks and Examples to appear with italic labels and an upright roman body. 

\begin{theorem} 

Sample theorem included for illustration. 

Numbers and parentheses, like equation $(3. 2)$, should be set 

in roman type. Note that words (as opposed to ' 'log-like'' 

functions) in displayed equations, such as 

$$ x~2 = Y~2 \sin z~2 \mbox{ for all } x $$ 

will appear in italic type in a theorem, though normally 

they should appear in roman. \end{theorem} 

This sample produces Theorem 4.1 below. 

Theorem 5.1. Sample theorem included for illustration. Numbers and parenthe- 
ses, like equation (3.2), should be set in roman type. Note that words (as opposed to 
"log-like" functions) in displayed equations, such as 

x 2 = Y 2 sin z 2 for all x 

will appear in italic type in a theorem, though normally they should appear in roman. 

Proofs are handled with the \begin{proof } \end{proof } environment. A "QED" 
box □ is created automatically by \end{proof }, but this should be preceded with a 
\qquad. 
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Named proofs, if used, must be done independently by the authors. SIAM style 
specifies that proofs which end with displayed equations should have the QED box 
two ems (\qquad) from the end of the equation on line with it horizontally. Below is 
an example of how this can be done: 

{\em Proof}. Proof of the previous theorem 



thus , 
$$ 

a~2 + b~2 = c~2 \qquad\endproof 
$$ 



6. Figures and tables. Figures and tables sometimes require special consid- 
eration. Tables in SIAM style are need to be set in eight point size by using the 
\f ootnotesize command inside the \begin{table} environment. Also, they should 
be designed so that they do not extend beyond the text margins. 

SIAM style requires that no figures or tables appear in the references section of 
the paper. I^TgX is notorious for making figure placement difficult, so it is important 
to pay particular attention to figure placement near the references in the text. All 
figures and tables should be referred to in the text. 

SIAM supports the use of epsf ig for including PostScript figures. All Post- 
Script figures should be sent in separate files. See the epsf ig documentation (avail- 
able via anonymous ftp from CTAN: ftp.shsu.edu) for more details on the use of this 
style option. It is a good idea to submit high-quality hardcopy of all PostScript 
figures just in case there is difficulty in the reproduction of the figure. Figures pro- 
duced by other non-TjrX methods should be included as high-quality hardcopy when 
the manuscript is submitted. 

PostScript figures that are sent should be generated with sufficient line thick- 
ness. Some past figures authors have sent had their line widths become very faint 
when SIAM set the papers using a high-quality 1200dpi printer. 

Hardcopy for non-PoSTSCRiPT figures should be included in the submission of 
the hardcopy of the manuscript. Space should be left in the {figure} command for 
the hardcopy to be inserted in production. 

7. Bibliography and BibTjrpC. If using BibTjtX, authors need not submit 
the .bib file for their papers. Merely submit the completed .bbl file, having used 
siam.bst as their bibliographic style file, siam.bst only works with BibTjrX version 
99i and later. The use of BibTgX and the preparation of a .bib file is described in 
greater detail in ||. 

If not using BibTjrX, SIAM bibliographic references follow the format of the fol- 
lowing examples: 

\bibitem{Ri} {\sc W. Riter}, 

{\em Title of a paper appearing in a book}, in The Book 
Title, E.~D. One, E.~D. Two, and A.~N. Othereditor, eds . , 
Publisher, Location, 1992, pp. "000—000. 

\bibitem{AuThl} {\sc A.~U. Thorone}, {\em Title of paper 
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with lower case letters}, SIAM J. Abbrev. Correctly, 2 
(1992), pp. -000—000. 

\bibitem{AlA2> {\sc A.~U. Thorone and A.~U. Thortwo}, {\em 
Title of paper appearing in book}, in Book Title: With All 
Initial Caps, Publisher, Location, 1992. 

\bibitem{AlA22} \sameauthor, % generates the 3 em rule 
{\em Title of Book{\rm :} Note Initial Caps and {\rm ROMAN 
TYPE} for Punctuation and Acronyms}, Publisher, 
Location, pp. "000— 000, 1992. 

\bibitem{AuTh3} {\sc A.~U. Thorthree}, {\em Title of paper 
that's not published yet}, SIAM. J. Abbrev. Correctly, to appear. 

Other types of references fall into the same general pattern. See the sample 
file or any SIAM journal for other examples. Authors must correctly format their 
bibliography to be considered as having used the macros correctly. An incorrectly 
formatted bibliography is not only time-consuming for SIAM to process but it is 
possible that errors may be introduced into it by keyboarders/copy editors. 

As an alternative to the above style of reference, an alphanumeric code may be 
used in place of the number (e.g., [AUTh90]). The same commands are used, but 
\bibitem takes an optional argument containing the desired alphanumeric code. 

Another alternative is no number, simply the authors' names and the year of 
publication following in parentheses. The rest of the format is identical. The macros 
do not support this alternative directly, but modifications to the macro definition are 
possible if this reference style is preferred. 

8. Conclusion. Many other style suggestions and tips could be given to help au- 
thors but are beyond the scope of this document. Simple mistakes can be avoided by 
increasing your familiarity with how IM^X functions. The books referred to through- 
out this document are also useful to the author who wants clear, beautiful typography 
with minimal mistakes. 

Appendix. The use of appendices. The \appendix command may be used 
before the final sections of a paper to designate them as appendices. Once \appendix 
is called, all subsequent sections will appear as 

Appendix A. Title of appendix. Each one will be sequentially lettered instead 
of numbered. Theorem-like environments, subsections, and equations will also have 
the section number changed to a letter. 

If there is only one appendix, however, the \Appendix (with a capital letter) 
should be used instead. This produces only the word Appendix in the section title, 
and does not add a letter. Equation numbers, theorem numbers and subsections of 
the appendix will have the letter "A" designating the section number. 

If you don't want to title your appendix, and just call it Appendix A. for 
example, use \appendix\section*{} and don't include anything in the title field. 
This works opposite to the way \section* usually works, by including the section 
number, but not using a title. 

Appendices should appear before the bibliography section, not after, and any 
acknowledgments should be placed after the appendices and before the bibliography. 
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PAUL DUGGAN+ AND VARIOUS A. U. THORS* 

Abstract. An example of SIAM IATJ5X macros is presented. Various aspects of composing 
manuscripts for SIAM's journal series are illustrated with actual examples from accepted manuscripts. 
SIAM's stylistic standards are adhered to throughout, and illustrated. 

Key words, sign-nonsingular matrix, LU-factorization, indicator polynomial 

AMS subject classifications. 15A15, 15A09, 15A23 

1. Introduction and examples. This paper presents a sample file for the use 
of SIAM's fATgX macro package. It illustrates the features of the macro package, 
using actual examples culled from various papers published in SIAM's journals. It is 
to be expected that this sample will provide examples of how to use the macros to 
generate standard elements of journal papers, e.g., theorems, definitions, or figures. 
This paper also serves as an example of SIAM's stylistic preferences for the formatting 
of such elements as bibliographic references, displayed equations, and equation arrays, 
among others. Some special circumstances are not dealt with in this sample file; for 
such information one should see the included documentation file. 

Note: This paper is not to be read in any form for content. The conglomeration of 
equations, lemmas, and other text elements were put together solely for typographic 
illustrative purposes and don't make any sense as lemmas, equations, etc. 

1.1. Sample text. Let S = [sy] (1 < i,j < n) be a (0, 1, — l)-matrix of order n. 
Then S is a sign-nonsingular matrix (SNS-matrix) provided that each real matrix with 
the same sign pattern as S is nonsingular. There has been considerable recent interest 
in constructing and characterizing SNS-matrices jl| , Q . There has also been interest 
in strong forms of sign-nonsingularity In this paper we give a new generalization 
of SNS-matrices and investigate some of their basic properties. 

Let S = [sij] be a (0, 1, — l)-matrix of order n and let C = [cy] be a real matrix 
of order n. The pair (S, C) is called a matrix pair of order n. Throughout, X = [xij] 
denotes a matrix of order n whose entries are algebraically independent indeterminates 
over the real field. Let S o X denote the Hadamard product (entrywise product) of 
S and X. We say that the pair (5*, C) is a sign-nonsingular matrix pair of order n, 
abbreviated SNS-matrix pair of order n, provided that the matrix 

A = SoX + C 

is nonsingular for all positive real values of the Xij . If C — O then the pair (S, O) is 
a SNS-matrix pair if and only if S is a SNS-matrix. If S = O then the pair (O, C) 
is a SNS-matrix pair if and only if C is nonsingular. Thus SNS-matrix pairs include 
both nonsingular matrices and sign-nonsingular matrices as special cases. 
The pairs (S, C) with 



S = 



1 




C = 



1 1 
1 1 
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are examples of SNS-matrix pairs. 

1.2. A remuneration list. In this paper we consider the evaluation of integrals 
of the following forms: 



(1.1) 



(1.2) 



J f(t) (j2EiB iXx (tyj dt, 



where Bi^^x is the ith B-spline of order k defined over the knots Xi,Xi+\, . . . ,Xi+k- 
We will consider B-splines normalized so that their integral is one. The splines may 
be of different orders and defined on different knot sequences x and y. Often the limits 
of integration will be the entire real line, — oo to +oo. Note that (1.1) is a special case 



of (1.2) where f(t) is a spline. 

There are five different methods for calculating ( jl.lD that will be considered: 

1. Use Gauss quadrature on each interval. 

2. Convert the integral to a linear combination of integrals of products of B- 
splines and provide a recurrence for integrating the product of a pair of B-splincs. 

3. Convert the sums of B-splines to piecewise Bezier format and integrate seg- 
ment by segment using the properties of the Bernstein polynomials. 

4. Express the product of a pair of B-splines as a linear combination of B- 
splines. Use this to reformulate the integrand as a linear combination of B-splines, 
and integrate term by term. 

5. Integrate by parts. 



Of these five, only methods 1 and 5 are suitable for calculating (1.2). The first four 
methods will be touched on and the last will be discussed at length. 

1.3. Some displayed equations and {eqnarray}s. By introducing the prod- 
uct topology on R mxm x R nxn with the induced inner product 

(1.3) ((AuBt), (A 2 , B 2 )) := (A ± ,A 2 ) + (B 1 ,B 2 ), 

we calculate the Frechet derivative of F as follows: 



F'(U, V)(H, K) = (R(U, V),HYV T + UY,K T - P(HY,V T + UT,K T )) 



(1.4) 



= {R(U,V),HZV T 
= (R(U,V)VX T ,H) 



^U 1 R{U,V),K A 



In the middle line of ( jl.4| ) we have used the fact that the range of R is always per- 
pendicular to the range of P. The gradient V-F of F, therefore, may be interpreted 
as the pair of matrices: 



(1.5) 



VF(C7, V) = (R(U, V)VE T , R(U, V) T UT.) e R mxm x R r > 
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Because of the product topology, we know 

(1.6) T (UiV) (0(m) x 0{n)) = TuO{m) x T v O(n), 

where F(ijy) (0(m) x 0(n)) stands for the tangent space to the manifold 0(m) x 0(n) 
at (U, V) G 0{m) x 0(n) and so on. The projection of VF(U, V) onto T {uy) (0(m) x 
0(n)), therefore, is the product of the projection of the first component of VF{U 1 V) 
onto TjjO(m) and the projection of the second component of VF(U, V) onto TvO(n). 
In particular, we claim that the projection g(U, V) of the gradient VF([7, V) onto 
T(ijy}(0(m) x 0(n)) is given by the pair of matrices: 



(1.7) 



g{Uj v) _^ R(u, v)v^ T u T uw T R(u, v) T L , 

R(U, V) T U^V T - VY?U T R{U, V) 



Thus, the vector field 



d(U, V) 



(1.8) ^ = - 9 (U,V) 

defines a steepest descent flow on the manifold 0(m) x 0(n) for the objective function 
F(U,V). 

2. Main results. Let (S, C) be a matrix pair of order n. The determinant 

det^oX + C*) 

is a polynomial in the indeterminates of X of degree at most n over the real field. We 
call this polynomial the indicator polynomial of the matrix pair (£>, C) because of the 
following proposition. 

Theorem 2.1. The matrix pair (S,C) is a SNS-matrix pair if and only if all 
the nonzero coefficients in its indicator polynomial have the same sign and there is at 
least one nonzero coefficient. 

Proof. Assume that (S, C) is a SNS-matrix pair. Clearly the indicator polynomial 
has a nonzero coefficient. Consider a monomial 



(2.1) h 



occurring in the indicator polynomial with a nonzero coefficient. By taking the Xij 
that occur in (2.1) large and all others small, we see that any monomial that occurs 
in the indicator polynomial with a nonzero coefficient can be made to dominate all 
others. Hence all the nonzero coefficients have the same sign. The converse is im- 
mediate. □ 

For SNS-matrix pairs (S, C) with C = O the indicator polynomial is a homo- 



geneous polynomial of degree n. In this case Theorem 2.1 is a standard fact about 
SNS-matrices. 

Lemma 2.2 (Stability). GivenT > 0, suppose that ||e(t)||i )2 < h q ~ 2 forO <t <T 
and q > 6. Then there exists a positive number B that depends on T and the exact 
solution tp only such that for all <t <T, 

(2-2) |lk(*)lli, a <^*- 3 / 2 + ||e(t)|| lia ). 
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The function B(T) can be chosen to be nondecreasing in time. 

Theorem 2.3. The maximum number of nonzero entries in a SNS-matrix S of 
order n equals 

n 2 + 3n - 2 



with equality if and only if there exist permutation matrices such that P\S\Q = T n 
where 



(2.3) 



T — 



1 1 
1 1 
1 








1 1 1 
1 1 1 
1 1 1 



1 1 1 

1 1 



We note for later use that each submatrix of T n of order n — 1 has all Is on its 
main diagonal. 

We now obtain a bound on the number of nonzero entries of S in a SNS-matrix 
pair (S, C) in terms of the degree of the indicator polynomial. We denote the strictly 
upper triangular (0,l)-matrix of order m with all Is above the main diagonal by U m . 
The all Is matrix of size m by p is denoted by J m . p . 

PROPOSITION 2.4 (Convolution theorem). Let 



Then 



a*u(t)= / a(t — r)u(r)rfr, £€(0, oo). 
Jo 



a * u(s) = a(s)u(s). 



Lemma 2.5. For s > 0, if 

/•OO 

/ e-^V 1 ^)^)^ <0, 
Jo 



then 



r°° 1 

/ e- 2sat v 2 {t)dt < —v 2 (0). 
Jo 2s o 



Proof. Applying integration by parts, we obtain 



r e- 2sot [v 2 (t)-v 2 (0)}dt= lim ( —e- 2aot v 2 (t)) + — [°° e~ 2sot v^ (t)v(t)dt 

Jo V 2s o / s Jo 

< _!_ / e - 2sot v^(t)v(t)dt < o. 

s o Jo 



Thus 



Jo 



e- 2sot v 2 {t)dt <v 



: (o)/ 

Jo 



e -2s t dt = w 2( ). 

2s 
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Corollary 2.6. Let E satisfy (5) -(6) and suppose E satisfies (7) and (8) with 
a general G. Let G = Vxl + Vp, p G Hq(Q). Suppose that Vp and V x 3? satisfy 
all the assumptions of Theorems 4.1 and 4.2, respectively. Ln addition suppose all the 
regularity assumptions of Theorems 4.1-4.2 are satisfied. Then for < t < T and 
< e < eo there exists a constant C — C(e, T) such that 

\\(E - E h )(t)\\ < Ch k+1 ~ € , 

where C also depends on the constants given in Theorems 4.1 and 4.2. 

Definition 2.7. Let S be an isolated invariant set with isolating neighborhood 
N . An index pair for S is a pair of compact sets {N\,Nq) with Nq C N± C N such 
that: 

(i) cl(Ni\N ) is an isolating neighborhood for S. 

(ii) iVj is positively invariant relative to N for i = 0, 1, i.e., given x € iVj and 
x ■ [0, t] C N, then x ■ [0, t] C N t . 

(iii) iVo is an exit set for N±, i.e. if x S N±, x ■ [0, oo) (t N±, then there is a 
T > such that x ■ [0, T] C JVi and x - T e N . 

2.1. Numerical experiments. We conducted numerical experiments in com- 
puting inexact Newton steps for discretizations of a modified Bratu problem, given 

by 

Aw + ce w + d^- = f in D, 

(2.4) dx 

w = on dD, 

where c and d are constants. The actual Bratu problem has d — and / = 0. It 
provides a simplified model of nonlinear diffusion phenomena, e.g., in combustion and 
semiconductors, and has been considered by Glowinski, Keller, and Rheinhardt 
as well as by a number of other investigators; see |jTT[| and the references therein. See 
also problem 3 by Glowinski and Keller and problem 7 by Mittelmann in the collection 
of nonlinear model problems assembled by More |ll| . The modified problem ( |2.4| ) has 
been used as a test problem for inexact Newton methods by Brown and Saad |7J . 
In our experiments, we took D = [0, 1] X [0, 1], / = 0, c = d — 10, and discretized 



( 2A ) using the usual second-order centered differences over a 100 x 100 mesh of equally 
spaced points in D. In GMRES(m), we took m = 10 and used fast Poisson right 
preconditioning as in the experiments in §2. The computing environment was as 
described in §2. All computing was done in double precision. 

In the first set of experiments, we allowed each method to run for 40 GMRES(m) 
iterations, starting with zero as the initial approximate solution, after which the limit 
of residual norm reduction had been reached. The results are shown in Fig. |2.1[ In 



Fig. 2.1, the top curve was produced by method FD1. The second curve from the top 
is actually a superposition of the curves produced by methods EHA2 and FD2; the 
two curves are visually indistinguishable. Similarly, the third curve from the top is 
a superposition of the curves produced by methods EHA4 and FD4, and the fourth 
curve from the top, which lies barely above the bottom curve, is a superposition of 
the curves produced by methods EHA6 and FD6. The bottom curve was produced 
by method A. 

In the second set of experiments, our purpose was to assess the relative amount of 
computational work required by the methods which use higher-order differencing to 
reach comparable levels of residual norm reduction. We compared pairs of methods 
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Fig. 2.1. Logio of the residual norm versus the number o/ GMRES(m) iterations for the finite 
difference methods. 

Table 2.1 

Statistics over 20 trials o/GMRES(m) iteration numbers, F -evaluations, and run times required 
to reduce the residual norm by a factor of e. For each method, the number o/ GMRES(m) iterations 
and F -evaluations was the same in every trial. 







Number of 


Number of 


Mean Run Time 


Standard 


Method 


e 


Iterations 


.F-Evaluations 


(Seconds) 


Deviation 


EHA2 


10 -io 


26 


32 


47.12 


.1048 


FD2 


io- 10 


26 


58 


53.79 


.1829 


EHA4 


io- 12 


30 


42 


56.76 


.1855 


FD4 


io- 12 


30 


132 


81.35 


.3730 


EHA6 


IO" 12 


30 


48 


58.56 


.1952 


FD6 


IO' 12 


30 


198 


100.6 


.3278 



EHA2 and FD2, EHA4 and FD4, and EHA6 and FD6 by observing in each of 20 
trials the number of GMRES(to) iterations, number of ^-evaluations, and run time 
required by each method to reduce the residual norm by a factor of e, where for 
each pair of methods e was chosen to be somewhat greater than the limiting ratio 
of final to initial residual norms obtainable by the methods. In these trials, the 
initial approximate solutions were obtained by generating random components as 
in the similar experiments in §2. We note that for every method, the numbers of 
GMRES(to) iterations and F'-evaluations required before termination did not vary at 
all over the 20 trials. The GMRES(m) iteration counts, numbers of J 71 - evaluations, 



and means and standard deviations of the run times are given in Table 2.1 



In our first set of experiments, we took c = d = 10 and used right preconditioning 
with a fast Poisson solver from FISHPACK which is very effective for these 
fairly small values of c and d. We first started each method with zero as the initial 
approximate solution and allowed it to run for 40 GMRES(m) iterations, after which 
the limit of residual norm reduction had been reached. Figure 2.2 shows plots of the 
logarithm of the Euclidean norm of the residual versus the number of GMRES(m) 
iterations for the three methods. We note that in Fig. 2/2 and in all other figures 
below, the plotted residual norms were not the values maintained by GMRES(m), 
but rather were computed as accurately as possible "from scratch." That is, at each 
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Fig. 2.2. Logio of the residual norm versus the number of GMRES(m) iterations for c = d = 10 
with fast Poisson preconditioning. Solid curve: Algorithm EHA; dotted curve: FDP method; dashed 
curve: FSP method. 



GMRES(m) iteration, the current approximate solution was formed and its product 
with the coefficient matrix was subtracted from the right-hand side, all in double 
precision. It was important to compute the residual norms in this way because the 
values maintained by GMRES(m) become increasingly untrustworthy as the limits 
of residual norm reduction are neared; see [|l7|]. It is seen in Fig. 2.2 that Algorithm 



EHA achieved the same ultimate level of residual norm reduction as the FDP method 
and required only a few more GMRES(m) iterations to do so. 

In our second set of experiments, we took c = d = 100 and carried out trials 
analogous to those in the first set above. No preconditioning was used in these ex- 
periments, both because we wanted to compare the methods without preconditioning 
and because the fast Poisson preconditioning used in the first set of experiments is not 
cost effective for these large values of c and d. We first allowed each method to run 
for 600 GMRES(m) iterations, starting with zero as the initial approximate solution, 
after which the limit of residual norm reduction had been reached. 

Acknowledgments. The author thanks the anonymous authors whose work 
largely constitutes this sample file. He also thanks the INFO-TeX mailing list for the 
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